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ABSTRACT 


Composite  strength  is  an  inverse  function  of  the  size  of  the  com¬ 
posite.  As  the  use  of  composites  expands  into  larger  applications, 
such  as  airplane  wings,  missile  components,  and  ship  superstructures, 
the  ability  to  accurately  predict  composite  performance  for  large 
applications  has  become  more  important.  The  composite  failure  pro¬ 
cess  is  sequential  and  initiates  with  early  breaking  of  the  weak  fibers. 
Concentration  of  breakage  sites  accumulates  and  leads  to  ultimately 
catastrophic  failure.  Prestressing  fibers  prior  to  solidification  of  the 
matrix  has  been  demonstrated  to  increase  the  reliability  of  the  com¬ 
posite  by  minimizing  the  spatial  concentration  of  the  breakage  sights. 

This  study  concentrates  on  quantifying  the  level  of  preload  and 
gauge  length  to  optimize  the  prestress  effect.  Computer  simulations  of 
graphite  bundle  tests  were  used  to  validate  the  data  analysis  method¬ 
ologies  applied  to  actual  AS-4  graphite  bundle  tests.  The  actual 
experimental  results  are  consistent  with  computer-simulated 
behaviors. 
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I.  INTRODUCTION 


In  the  1980s,  the  capability  of  military  equipment  has  increased 
exponentially.  Modem  aircraft  fly  further  and  carry  greater  payloads, 
missiles  have  larger  warheads  and  greater  ranges,  tanks  have  bigger 
guns  and  tougher  armor,  and  ships  have  more  capable  but  lighter 
weapon  systems  along  with  improved  armor  plating.  Engineering 
improvements  in  the  areas  of  electronics  and  composite  materials 
have  been  on  the  leading  edge  of  this  technological  revolution. 

As  reported  in  the  fourteenth  edition  of  Ships  and  Aircraft  of  the 
U.  S.  Fleet,  the  AV-8B  Harrier  II  VSTOL  attack  aircraft  “provides  twice 
the  payload  of  the  AV-8A  with  up  to  9200  lbs  of  external  stores."  The 
primary  difference  between  the  AV-8A  and  the  AV-8B  is  the  extensive 
use  of  composite  materials  in  the  wing  and  fuselage.  This  is  just  one 
example  of  the  use  of  composite  materials  in  existing  military  applica¬ 
tions.  New  applications  include  the  VS-22  Osprey  tilt  rotor  aircraft, 
whose  wing  and  fuselage  will  be  entirely  made  of  graphite  epoxy  com¬ 
posite  material.  Emerging  “stealth”  technology  and  its  application  to 
aircraft  is  a  further  new  application  of  composites. 

Composites  offer  designers  many  attractive  features.  These 
include  very  high  strength-to-weight  ratios,  high  stiffness-to-weight 
ratios,  wear  resistance,  controllable  heat  expansion  coefficients,  cor¬ 
rosion  resistance,  fatigue  resistance,  low  price,  and  ease  of  manufac¬ 
turing  and  repairability.  The  field  is  expanding  rapidly  and  research  is 


bringing  to  the  designer  stronger,  tougher,  and  stiffer  materials  every 
day. 

A  hurdle  to  be  overcome  for  composites,  however,  is  the  general 
lack  of  understanding  of  the  physics  and  mechanisms  at  work  in  the 
composite  which  give  these  materials  their  unique  properties.  This 
lack  of  understanding  has  restrained  designers’  ability  to  accurately 
predict  composite  strength.  To  compensate  for  this,  large  factors  of 
safety  are  being  used  in  the  design  of  composite  structures,  which  in 
turn  adds  weight  and  size  and  causes  other  related  problems  that 
reduce  design  efficiency.  As  the  application  of  composite  materials  has 
expanded  to  larger  and  larger  components,  such  as  entire  fuselages, 
rocket  motor  cases,  and  gun  turrets,  this  problems  has  become 
magnified. 

Composite  materials  are  two  or  more  dissimilar  materials  that 
when  properly  combined  form  a  new  material  whose  properties 
exceed  either  of  the  individual  constituents.  For  the  designer,  com¬ 
posite  materials  provide  an  ability  not  only  to  optimize  his  design  but 
also  to  optimize  the  material  to  be  used  for  the  design.  Material  prop¬ 
erties  of  conventional  materials  that  can  be  improved  upon  by  using 
composites  include  stiffness,  corrosion  resistance,  wear  resistance, 
strength-to-weight  ratio,  and  radar  wave  absorbability  and  reflectivity. 
Composites  in  use  today  include  fiberglass,  WEST  (Wood  Epoxy 
Saturation  Technique),  graphite-epoxy,  Aramid-epoxy,  and  Boron- 
epoxy.  Composite  materials  also  include  graphite-aluminum  and  hybrid 
composites  of  more  than  two  constituents. 
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Commonly  accepted  types  of  composites  are  the  fibrous  compos¬ 
ites,  which  consist  of  fibers  in  a  matrix;  laminated  composites,  which 
consist  of  layers  of  various  materials;  and  particulate  composites,  which 
are  composed  of  particles  in  a  matrix.  The  composite  is  fabricated  by 
combing  the  reinforcement  material  into  a  generally  ductile  matrix 
such  as  epoxy.  This  lay-up  is  called  a  lamina.  Several  layers  of  lamina 
may  be  combined  in  prescribed  geometrical  orientations  to  create  a 
laminate.  The  designer  specifies  his  materials  and  then  designs  the 
laminate  to  meet  his  requirements. 

Until  the  early  1980s  material  properties  of  composites  were 
predicted  by  a  simple  principle  called  the  “rule  of  mixtures."  This 
approximation  assumed  that  each  constituent  of  the  composite  con¬ 
tributed  to  the  overall  composite  property  in  direct  proportion  to  its 
property  weighted  by  its  volume  percentage  in  the  composite.  This 
linear  weighting  method  was  satisfactory  for  composite  physical  prop¬ 
erties.  which  are  averaged  combinations  of  the  constituent  properties 
such  as  density  and  stiffness. 

However,  the  "rule  of  mixture”  is  very  inaccurate  for  those  physi¬ 
cal  properties  which  are  governed  by  extreme  values  such  as  perme¬ 
ation  and  strength.  For  such  properties,  considerable  testing  of  a 
specific  composite  is  needed  to  achieve  representative  material  prop¬ 
erty  values  suitable  for  engineering  design  and  reliability  assurance. 
Not  only  are  large  amounts  of  data  required  for  each  specific 
composite  but  the  designer's  material  selections  become  limited  to 
only  those  composites  which  have  been  tested.  This  characterization 


3 


methodology  could  not  keep  up  with  the  rapid  pace  of  the 
introduction  of  new  materials. 

Clearly,  a  characterization  method  was  needed  to  be  able  to 
predict  composite  material  properties,  particularly  strength,  if  the  full 
advantage  of  composites  was  to  be  obtained.  Much  work  has  been 
completed  in  the  area  by  researchers  such  as  Phoenix,  Rosen,  Harlow, 
and  Wu.  The  most  current  papers  on  reliability  and  strength  predic¬ 
tions  of  composites  may  be  found  in  References  1,  2,  3,  and  4.  The 
most  significant  part  of  their  efforts  is  that  they  have  identified  that 
the  strength  of  composites  is  a  probabilistic  manifestation  of  the  sta¬ 
tistical  strength  of  the  fiber  strength.  The  statistics  of  the  fiber 
strength  govern  the  number  of  failure  sites  in  a  composite.  The  proba¬ 
bilistic  spatial  distribution  of  the  failure  sites  together  with  the  failure 
mechanism  govern  the  strength  of  a  specific  composite. 

Different  probabilistic  models  have  evolved  to  capture  the  essence 
of  the  failure  mechanisms.  The  simplest  model  is  the  “weakest  link 
model,”  which  envisions  the  entire  composite  as  a  long  thin  fiber;  the 
entire  chain  (the  composite)  fails  when  the  weakest  link  of  the  fiber 
fails.  This  model  does  not  account  for  the  redundancy  provided  by  the 
matrix  and  as  a  result  grossly  underestimates  the  strength  of  the  com¬ 
posite.  The  “equal  load-sharing  model”  accounts  for  the  redundancy 
provided  by  the  matrix  through  definition  of  an  effective  length,  an 
isolation  parameter  for  the  initial  breaks.  However,  this  model  does 
not  account  for  the  stress  concentration  around  the  breaks  and,  as  a 
result,  over-predicts  the  composite  strength.  Finally,  the  “local 
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load-sharing  model”  accounts  for  both  the  matrix  redundancy  and  the 
local  stress  concentration.  These  refinements  lead  to  the  most 
realistic  predictions  to  date. 

There  are  two  phenomena  relating  to  the  statistical  strength  of 
composites.  The  first  is  that  the  larger  a  composite  becomes,  the 
weaker  it  gets.  The  second  phenomenon  is  that  the  larger  a  compos¬ 
ite.  the  less  strength  scatter  among  such  large  structures.  This  is  due 
to  the  the  lower  weak  tail  of  the  strength  distributions  of  the  con¬ 
stituent  fibers.  This  effect  has  been  demonstrated  in  the  laboratory 
and  in  the  field. 

An  extension  of  this  theory  then  becomes  that  if  the  fibers  are 
preloaded  to  break  the  weak  sections  of  the  fiber  prior  to  their  lay-up 
into  the  matrix,  then  the  lower  tail  of  the  composite  strength 
distribution  is  reduced  and.  in  fact,  becomes  bounded  by  a  value.  This 
process  is  called  prestressing.  For  a  prescribed  value,  the  strength  and 
reliability  of  the  composite  may  then  be  theoretically  calculated  for  a 
given  fiber  and  the  magnitude  of  preload  based  on  the  modified  fiber 
distribution.  The  variables  to  be  determined  to  prestress  a  composite 
include  the  magnitude  of  the  preload,  the  gauge  length  to  which  the 
load  is  applied,  and  the  modified  fiber  parameters. 

The  ability  to  accurately  recover  the  fiber  strength  parameters  for 
small  (less  than  2.5  cm)  gauge  lengths  is  very  important  if  the  pre¬ 
stress  effect  is  to  be  optimized.  The  bundle  testing  done  in  the  labora¬ 
tory  yields  data  that  includes  not  only  the  fiber  data  but  also  machine 
compliance,  slack,  and  friction.  For  short  gauge  lengths,  the 
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contributions  of  slack  and  friction  are  small  but  compliance  becomes  j 

relatively  large.  The  portion  of  displacement  contributed  to  by  the  M 

experimental  system  compliance  in  a  bundle  test  needs  to  be  removed 
from  the  data  if  the  true  fiber  parameters  are  to  be  obtained. 

This  study  entailed  the  mathematical  simulation  modeling  of  ■ 

strength  prediction  for  a  bundle  of  fibers  with  no  matrix.  This  simula¬ 
tion  was  run  on  an  IBM  Personal  Computer  using  Microsoft  Fortran 
4.01  for  source  code  and  Lotus  1-2-3  for  graphing.  When  the  Simula-  fl 

tion  program  had  been  validated,  a  zero  gauge  length  bundle  test  was 
completed  to  measure  the  INSTRON  fiber  bundle  testing  machine 
compliance.  A  set  of  procedures  was  written  to  optimize  a  curve  fit  of  1 

the  compliance  data  to  provide  an  explicit  equation  which  models  the 
INSTRON  system  compliance. 

When  this  procedure  had  been  validated,  bundle  tests  of  varying  f 

lengths  were  tested  at  the  Mechanics  of  Materials  for  Composites  Lab¬ 
oratory  at  the  Naval  Postgraduate  School.  The  displacement  con¬ 
tributed  to  the  test  data  by  the  INSTRON  system  compliance  was  j 

removed  from  the  total  displacement.  Under  another  thesis  effort,  this 
rectified  data  was  analyzed  in  order  to  extract  the  fiber  statistical 
strength  distribution  parameters.  < 

The  purpose  of  this  investigation  was  to  validate  the  experimental 
procedures  and  data  analysis  methods  with  the  use  of  computer  simu¬ 
lated  data.  Data  generated  from  Monte  Carlo  simulation  using  known  j 

parameters  is  then  utilized  to  assist  in  experimental  design  and  the 
associated  prediction  theory.  Experimental  design  and  data 
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interpretation  methodology  provide  a  proven  foundation  toward  iden¬ 
tification  of  an  optimum  gauge  length  for  prestressing  in  order  to 
guarantee  the  reliability  of  composite  structures. 


n.  BACKGROUND 


The  tensile  failure  process  of  composite  materials  is  a  complex 
sequential  combination  of  many  different  processes.  These  processes 
are  controlled  by  the  statistical  strength  of  the  fiber  and  the  relative 
spatial  clustering  of  the  broken  fiber  sites  within  a  composite.  Given  a 
probabilistic  model  which  is  capable  of  characterizing  the  sequential 
failure  mechanisms,  a  formula  can  be  derived  to  predict  the 
probabilistic  strength  of  the  composite  in  terms  of  the  statistical 
strength  of  the  fiber.  The  parameters  for  such  a  probabilistic  models 
can  then  be  expressed  in  terms  of  other  common  measures  of  central 
tendencies  (e.g..  mean  composite  strength)  and  dispersions  (e.g., 
strength  scatter). 

A.  COMPOSITE  STRENGTH  PREDICTION 

The  prediction  of  the  probabilistic  strength  of  a  composite  as  a 
function  of  fiber  statistics  was  first  introduced  by  B.  W.  Rosen  in 
Reference  5,  in  the  early  sixties.  He  refined  the  simple  weakest  link 
model  (which  is  essentially  a  series  model)  by  accounting  for  the 
redundancy  introduced  by  the  presence  of  a  matrix  binder.  The  com¬ 
posite  geometry  he  examined  is  shown  below: 

•  a  single  layer  of  fibers 

•  evenly  distributed  throughout  the  matrix 

•  no  fibers  are  in  contact  with  other  fibers 

•  the  load  is  purely  axial  and  in  tension. 
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The  strength  of  a  single  lamina  of  composite  material  is  con¬ 
trolled  by  four  factors: 

•  The  strength  and  modulus  distributions  of  the  fiber. 

•  The  strength  and  modulus  distributions  of  the  matrix. 

•  The  interaction  between  the  fiber  and  the  matrix  with  respect  to 
load  sharing  and  stress  concentration. 

•  The  location  distribution  of  relative  strengths  of  the  fibers  with 
respect  to  the  location  distribution  of  relative  strengths  of  the 
matrix. 

When  a  composite  is  pulled  in  tension,  because  the  fiber  generally 
has  a  much  higher  modulus  as  compared  to  the  matrix  (e.g.,  graphite 
108  Psi  versus  epoxy  104  psi),  the  fibers  will  carry  most  of  the  load. 
The  load-carrying  capability  of  the  matrix  can  be  neglected. 

As  the  load  is  increased,  each  fiber  will  carry  an  equal  amount  of 
load  that  is  equal  to  the  applied  load  divided  by  the  number  of  fibers. 

Pf  =  Pt/n  (1) 

The  load  on  each  fiber  will  be  equal,  but  because  the  diameter  of 
the  fibers  is  not  constant,  the  stress  carried  at  each  point  in  the  fiber 
will  be  different.  The  stress  will  vary  inversely  to  the  diameter  of  the 
fiber.  Some  of  the  variations  are  regions  of  gradual  thinning,  others  are 
indentations  or  nicks.  These  areas  of  reduced  thickness  or  indenta¬ 
tions  or  nicks  are  called  “flaws. "  These  flaws  reduce  the  load-carrying 
capability  of  the  fiber  by  reducing  its  thickness  and  introducing  stress 
concentration  areas  in  the  fiber. 
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A  great  deal  of  study  has  been  done  in  identifying  and  quantifying 
the  effect  of  these  flaws  on  a  microscopic  scale.  The  work  of  Phoenix 
and  Harlow  proposed  that  on  a  macroscopic  scale,  if  the  flaws  are  dis¬ 
tributed  randomly  with  respect  to  location  on  the  fibers  and  normally 
with  respect  to  the  severity  of  the  flaw,  then  the  strength  of  a  fiber 
can  be  modeled  as  a  long  chain  of  segments.  The  strength  of  the  seg¬ 
ments  may  be  modeled  using  the  Weibull  distribution.  The  model 
selection  is  expanded  in  Appendix  A.  The  Weibull  distribution  is  an 
appropriate  model  for  the  probability  of  failure  of  a  chain  with  many 
links. 


F(x)  =  1  -  exp  l-(x/px)ot]  (2) 

Where  x  =  is  the  load  applied  to  the  fiber  (the  chain) 

a  =  is  the  shape  parameter  characterizing  the  variability 

px  =  is  the  scale  parameter  characterizing  central  location 
of  the  failure  strengths. 

The  Weibull  distribution  for  strength  predicts  the  “size  effect." 
That  is,  if  the  parameters  are  known  for  one  given  length,  the 
parameters  for  any  other  lengths  of  the  same  material  can  be  com¬ 
puted.  This  is  also  expanded  in  Appendix  A. 

As  the  load  increases  on  the  material,  the  stress  on  each  segment 
of  the  fiber  will  increase  as  well.  When  the  load  reaches  a  point  where 
the  resulting  stress  is  equal  to  the  ultimate  strength  of  a  segment,  that 
segment  of  fiber  will  break.  This  break  will  release  a  certain  amount  of 
strain  energy  into  the  composite.  It  will  also  cause  an  immediate  local 
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load  concentration  that  will  be  felt  by  the  adjacent  and  surrounding 
fibers. 

The  matrix,  via  its  shear  stress  carrying  capability,  will  attempt  to 
absorb  the  released  strain  energy  and  distribute  the  local  load  concen¬ 
tration  to  the  surrounding  fibers.  The  actual  values  of  the  stress  con¬ 
centration  factors  caused  by  the  fiber  break  are  not  known,  but  shear 
leg  analysis  yields  reasonable  approximations. 

Once  one  fiber  segment  breaks,  many  events  can  occur.  The  first 
is  that  the  matrix  is  strong  enough  to  absorb  the  released  strain 
energy  and  the  surrounding  fiber  elements  are  capable  of  carrying  the 
local  load  concentration.  If  this  occurs,  the  material  is  stable  and  the 
load  on  each  fiber  will  increase  to 

Pf  «  Pt/(1  -  n)  (3) 

Where  n  represents  the  number  of  broken  fibers. 

Another  possibility  is  that  the  matrix  can  not  absorb  the  released 
strain  energy.  In  this  case  the  matrix  will  split  as  it  fails  in  shear.  The 
material  will  continue  to  tear  until  its  radius  is  such  that  the  stress 
concentration  decreases  to  a  point  below  the  ultimate  strength  of  a 
section  of  the  matrix.  The  fibers  in  this  region  will  then  be  subjected 
to  additional  stress  as  the  matrix  pulls  away.  If  the  segments  of  the 
matrix  are  strong  enough  to  carry  the  additional  load,  the  material  will 
then  again  become  stable. 

A  further  possibility  is  that  the  adjacent  fiber  segment  may  not  be 
able  to  withstand  the  stress  caused  by  the  first  fiber  failure.  In  this 


case,  that  fiber  segment  will  break  and  cause  additional  strain  energy 
to  be  released.  This  reaction  is  called  auto  clustering.  The  process  will 
continue  until  it  is  interrupted  by  a  fiber  segment  that  is  strong 
enough  to  carry  the  stress  or  the  entire  piece  fails  catastrophically. 

Some  observations  of  the  description  above:  The  composite  may 
fail  due  to  either  the  local  strength  of  the  matrix  or  the  local  strength 
of  the  fiber  or  a  combination  of  both.  The  likely  initiation  of  the  failure 
process  will  occur  with  a  fiber  break  because  the  strength  scatter  or 
variance  of  the  fiber  is  typically  greater  than  that  of  the  matrix  giving 
rise  to  nucleation  sites  at  already  high  stress  regions.  The  chain  of 
events  that  occur  after  that  first  break  is  dependent  on  the  strength 
statistics  of  the  adjoining  matrix  and  fibers.  If  there  are  many  strong 
segments  surrounding  a  weak  segment  that  fails,  then  additional 
material  will  not  fail.  If,  however,  this  is  not  the  case,  then  further 
damage  occurs  which  may  lead  to  final  unstable  catastrophic  failure. 

From  the  previously  described  failure  process.  It  can  be  observed 
that  the  strength  of  the  composite  is  dependent  on  the  nucleating 
weak  segments  and  the  probability  of  there  being  another  weak  seg¬ 
ment  nearby.  Therefore,  the  shape  factor  of  the  fiber  strength  distri¬ 
bution  is  the  critical  element  that  controls  composite  strength  and 
reliability.  Because  of  the  mechanism  of  failure,  as  a  composite 
becomes  larger,  the  number  of  the  population  of  fiber  segments 
increases  and  the  likelihood  of  the  existence  of  very  weak  segments 
becomes  larger.  This  causes  the  strength  of  the  entire  composite  to 
decrease. 
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Another  observation  is  that  because  the  strength  of  the  composite 
is  controlled  by  the  weak  fiber  segments,  the  strength  may  not  be 
solved  for  explicitly  because  the  actual  value  of  the  weak  segments  is 
probabilistic.  The  strength  of  adjoining  segments  is  probabilistic,  as 
well. 

Prestressing  solves  two  critical  problems.  First,  it  sets  the  lower 
limit  of  the  strength  of  the  segments  by  eliminating  the  weak  seg¬ 
ments  by  causing  them  to  break  prior  to  the  application  of  matrix. 
Second,  prestressing  causes  the  weak  segments  to  be  broken  without 
local  load  stress  concentrations  affecting  the  adjacent  fibers  and  with¬ 
out  the  localized  dynamic  strain  energy  release  because  the  weak  seg¬ 
ments  are  broken  before  the  matrix  is  solidified.  This  results  in  auto 
clustering  being  minimized. 

The  critical  implementation  problem  is  then  to  find  the  “segment 
length”  to  ensure  that  the  majority  of  the  weak  segments  are  elimi¬ 
nated  and  then  determine  a  prestress  level.  The  prestress  level  may 
be  optimized  to  a  level  where  the  advantages  of  breaking  the  weak 
segments  is  balanced  by  the  introduction  of  clusters  to  the  matrix. 

B.  PRESTRESSING 

G.  J.  Mills  and  associates,  under  an  United  States  Air  Force-spon¬ 
sored  project,  conducted  prestress  tests  on  Boron  fibers  in  the  early 
1970s.  In  his  reports,  documented  in  References  5,  6,  and  7,  he 
detailed  how  his  work  demonstrated  the  feasibility  of  the  concept  and 
provided  data  from  early  experiments  that  supported  his  expectations. 
His  paper  did  not  attempt  to  model  or  predict  the  effect  of 


prestressing.  He  did  note  the  effect  of  varying  prestress  gauge  length 
and  that  as  the  gauge  length  decreased,  the  strength  of  the  fiber 
increased,  but  he  did  not  attempt  to  quantify  an  optimum  length.  He 
described  the  strength  increase  of  the  composite  purely  in  terms  of 
removal  of  weak  segment  and  reducing  the  tail-end  scatter  of  the 
fiber.  He  did  not  attempt  to  predict  the  performance  of  the  composite 
once  the  matrix  had  been  set.  He  did  report  that  the  strength  did 
increase. 

In  1986,  Lt.  David  Bell  and  Professor  E.  Wu  of  the  Naval  Postgrad¬ 
uate  School  conducted  prestress  tests  on  AS-4  graphite/epoxy.  The 
results  are  published  in  Reference  9.  The  report’s  experimental  data 
confirmed  the  phenomenon  of  strength  improvement  after  pre¬ 
stressing.  This  was  in  concurrence  with  Mills’  laboratory  data.  Both 
papers  accounted  for  the  improvement  by  theorizing  that  the  process 
of  prestressing  improved  the  strength  of  a  composite  because  it 
reduced  the  strength  scatter  of  the  composite.  The  work  done  by  Bell 
and  Wu  was  only  done  with  one  gauge  length,  10",  and  even  with  the 
problems  encountered  in  handling  the  samples,  the  material  did  show 
a  significant  reduction  in  scatter  and  therefore  an  increase  in  the 
usable  strength  reliability. 

Implementation  of  these  findings  to  a  practical  application 
requires  a  model  that  can  adequately  predict  the  effect  of  prestressing 
as  a  function  of  prestress  level,  gauge  length,  and  fiber  strength 
parameters. 
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C.  TESTING  METHODS 

There  are  two  methods  in  which  fiber  strength  statistics  can  be  I 

measured.  One  method  is  to  prepare  and  perform  a  breaking  strength 
test  on  one  individual  fibers.  The  second  method  is  to  form  many 
fibers  into  a  bundle  of  parallel  fibers  and  then  perform  a  controlled  j 

strain  rate  failure  test. 

The  first  method  is  very  time  consuming  and  error  prone.  The 
fibers  are  very  thin  (in  the  order  of  10  pm)  and  handling  requires  j 

delicate  treatment.  Laboratory  conditions  must  be  carefully  monitored 
as  the  small  fibers  are  susceptible  to  damage  caused  by  temperature, 
wind,  abrasion,  and  other  disturbances.  In  addition,  each  test  only  I 

provides  one  single  point  of  data.  If  failure  probability  data  is  required 
down  to  .001  failure  rate,  then  at  least  10,000  samples  should  be 
tested.  I 

Obviously,  the  second  method  is  the  preferred  method.  The  bun¬ 
dle  test,  performed  in  our  case  on  bundles  of  1,000  filaments,  will 
yield  1,000  data  points  with  one  test.  The  handling  problems,  while  I 

still  present,  are  on  a  much  more  manageable  level.  The  data  analysis 
of  the  load  vs.  displacement  data  will  yield  the  strength  distribution 
parameters  for  the  fiber  that  is  equal  to  1,000  individual  filament  i 

tests. 

It  is  important  then  that  we  measure  the  true  load  and  displace¬ 
ment  for  just  the  fibers.  The  displacement  that  the  testing  machine 
measures  is  the  sum  of  the  fiber  displacement,  machine  compliance, 
electronic  noise,  slack  in  the  individual  fibers,  and  friction  between 
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the  fibers.  Figure  1  demonstrates  that  test  data  is  the  combination  of 
all  of  these  factors. 


Figure  1 

Typical  Bundle  Test  Data  Graph 


The  compliance  is  due  to  the  mechanical  construction  of  the 
testing  machine.  The  machine  operates  with  gears,  clutches,  and 
motors.  Each  of  these  has  an  initial  start-up  slip  or  backlash.  The 
machine  also  has  grips  which  hold  the  samples  and  the  bases  which 
will  elongate  under  load.  They  are  very  stiff  and  the  elongation  is  small 
when  compared  to  the  bundles,  except  when  the  bundle  gauge  length 
is  short.  This  is  illustrated  in  Figure  2,  where  the  displacement  is 
shown  as  percent  strain.  The  compliance  for  the  long  sample  has  a 


greater  slope  than  the  short  sample  because  the  slope  of  the  compli¬ 
ance  in  displacement  multiplied  is  by  the  gauge  length  for  percent 
strain. 


%  Strain  3  %  Strain  3 


GL  =  Long  (>  25  cm)  GL  =  Short  (<  2.5  cm) 

Figure  2 

Effect  of  Compliance  on  Long-Gauge  Length 
vs.  Short-Gauge  Length 

The  compliance  of  the  machine  can  be  measured  if  the  there  is  no 
contribution  to  the  displacement  from  the  fibers.  This  is  achieved 
using  a  zero  gauge  length  sample.  A  sample  must  be  utilized  to  include 
the  sample  grip  mechanism  in  the  compliance.  A  curve-fitting 
approximation  is  then  applied  to  the  data  to  create  a  continuous  func¬ 
tion  of  d  in  terms  of  P.  This  function  may  then  be  subtracted  from  the 
test  data  to  eliminate  the  effect  of  compliance.  This  procedure  and 
testing  is  expanded  in  Appendix  C. 
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The  effect  of  noise  is  sensed  in  the  load  cell.  The  sensitivity  of  the 
load  cell  is  1  out  of  10,000.  This  is  due  to  the  limit  of  the  analog-to- 
digital  converter,  which  can  only  manage  four  bits  of  data.  The  amount 
of  noise  in  the  machine  sensed  by  the  load  cell  is  very  small  when 
compared  to  the  sensitivity  of  the  load  cell.  As  the  load  increase  and 
the  output  changes  in  order  of  magnitude,  the  noise  is  not  even 
noticeable.  The  noise  in  the  machine  is  considered  to  be  random 
electrical  white  noise  due  to  interference  in  the  data  transmission 
lines.  Observation  of  the  data  indicates  that  the  noise  is  not  very 
significant. 

A  simulation  program  was  written  during  this  study  that  allowed 
the  user  to  input  random  noise  up  to  any  level.  This  will  allow  simula¬ 
tions  to  be  run  and  investigate  the  effect  of  noise  on  the  fiber  strength 
parameters  recovered  from  the  simulated  test. 

The  slack  is  due  to  the  limitation  of  the  sample  preparation  that 
all  the  fibers  do  not  have  exactly  the  same  gauge  length.  This  problem 
is  found  to  be  greater  with  the  longer  samples  as  compared  to  the 
short  samples.  The  slack  can  not  be  eliminated  by  current  sample 
preparation  techniques,  but  extreme  care  by  the  technician  in 
preparing  the  sample  can  minimize  the  effect. 

Lt.  Joseph  Schmidt  is  examining  the  effects  of  slack  on  the  calcu¬ 
lated  fiber  parameters  and  expects  to  publish  the  results  of  his  study 
in  September  1988,  from  the  Naval  Postgraduate  School,  Monterey. 

The  effect  of  friction  in  the  bundle  test  cannot  not  be  definitively 
characterized.  The  presence  of  friction  is  manifested  by  large  vertical 
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drops  In  the  load  displacement  curve,  perhaps  due  to  the  tangling  of 
broken  fibers  with  unbroken  fibers.  The  friction  also  will  cause  a  fibers 
to  break  more  than  once.  The  effects  of  friction  do  increase  with 
longer  gauge  lengths  and  are  reduced  with  shorter  gauge  lengths. 

With  these  factors  that  affect  the  data  of  the  testing  machine 
accounted  for,  the  true  fiber  strength  statistics  may  be  recovered  from 
the  bundle  tests.  These  statistics  may  then  by  used  to  optimize  the 
prestress  effect  by  calculating  the  effect  prestressing  will  have  in 
altering  the  fiber  strength  distribution. 
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III.  PRESTRESSING  MECHANISM 


The  prestressing  process  Is  to  preload  the  composite  fibers  prior 
to  solidification  of  the  matrix  to  cause  the  weak  segments  of  the  fibers 
to  break.  The  load  that  was  carried  by  the  broken  filaments  is  then 
shared  equally  among  the  remaining  fibers.  Equal  load  sharing  is  in 
effect  because  the  pre-solidified  matrix  cannot  support  a  shear  load 
and  cause  stress  concentration  cells  to  form.  Upon  solidification  of  the 
matrix,  however,  local  load  sharing  is  in  effect,  thereby  increasing  the 
local  redundancy.  Because  the  weak  segments  of  the  fibers  were 
broken  prior  to  solidification  of  the  matrix,  the  local  load  concentra¬ 
tion  stress  cells  are  not  as  severe  as  if  the  fibers  had  broken  while  in 
the  matrix.  The  reduction  in  the  severity  of  the  stress  concentration 
cells  provides  the  reliability  enhancement  of  the  composite  and  is  the 
motivation  for  prestressing. 

The  parameters  that  affect  the  prestress  effect  are  gauge  length, 
prestress  magnitude,  and  the  fiber  strength  distribution  statistics.  The 
optimum  gauge  length  to  prestress  the  fibers  would  be  the  ineffective 
length.  The  ineffective  length,  as  described  by  Rosen  in  Reference  5, 
is  a  function  of  both  the  tensile  strength  of  the  fiber  and  of  the  shear 
strength  of  the  matrix.  This  makes  the  ineffective  length  unique  for 
each  composite,  and  because  both  strengths  are  random  variables,  the 
ineffective  length  for  a  particular  location  in  the  composite  is  random 
as  well. 
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As  the  gauge  length  increases,  the  likelihood  of  breaking  all  of  the 
weak  segments  in  the  fiber  decreases.  This  is  because  using  the  equiv¬ 
alent  load  sharing  model,  it  is  only  possible  to  cause  one  break  in  each 
fiber  over  the  gauge  length.  The  break  occurs  at  the  weakest  flaw, 
which  may  occur  before  reaching  the  desired  prestress  level.  As  illus¬ 
trated  in  Figure  3,  this  allows  other  flaws  weaker  than  the  prestress 
level  to  remain  unbroken. 


p  INEFFECTIVE 
L  LENGTH 


r)  PRESTRESS 
0L '  GAUGE  LENGTH 


GL 


Flaw 

Broken 


Flaw 

Not  Broken 


M  =  1  M  =  3  M  =  6 

Figure  3 

Effect  of  Preloading  for  Various  Gauge  Lengths 


| 


21 


When  the  prestressed  fiber  is  placed  into  a  matrix,  the  local  load- 
sharing  model  becomes  applicable.  With  this  model,  the  weak  flaws  in 
the  composite  will  break  under  load  and  generate  the  larger  stress 
concentration  cells  which  prestressing  is  designed  to  minimize.  The 
mathematical  modeling  of  solving  for  the  probability  density  function 
of  remaining  weak  flaws  given  gauge  lengths  greater  than  the  ineffec¬ 
tive  length  and  prestress  levels  is  very  complex.  This  problem  is  being 
treated  in  a  separate  project. 

The  preload  level  has  two  effects  on  the  prestress  effect.  For 
higher  prestress  levels,  more  fiber  breaks  are  introduced  and  more 
stress  concentration  cells  will  be  created  in  the  the  composite  once 
the  matrix  has  solidified.  The  severity  of  the  concentration  cells  is  a 
function  of  the  ineffective  length.  The  shorter  the  ineffective  length  is, 
the  smaller  the  stress  concentration  cell  is  because  of  less  local  load 
sharing.  Second,  the  higher  the  preload  level,  the  more  difficult  it 
becomes  to  grip  or  maintain  the  desired  gauge  length.  The  fibers  by 
their  nature  have  a  small  diameter,  on  the  order  of  5  microns  for 
graphite,  and  are  very  susceptible  to  handling  damage. 

The  underlying  fiber  statistics  also  have  an  effect  on  the  prestress 
effect.  The  fiber  strength  variability  (the  Weibull  model  shape  parame¬ 
ter  a),  has  the  greatest  influence  on  the  prestress  effect.  As  seen  in 
Figure  4  ,  if  the  shape  parameter  is  very  high  (a  >  20  ),  i.e.,  the  fiber 
has  low  variability,  then  the  prestress  effect  will  not  be  very  helpful. 
This  is  because  the  range  of  the  distribution  of  the  fiber  strengths  is 
very  small  and  the  weak  tail  that  prestressing  eliminates  or  reduces  is 
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Inherently  not  present.  Fibers  of  this  strength  uniformity,  however, 
are  not  available  at  this  time.  If  the  shape  parameter  is  small,  then  the 
scatter  is  much  greater.  With  this  scatter,  the  weak  tail  may  be  altered 
or  removed  by  prestressing. 


Figure  4 

Strength  Distribution  for  Fibers  (Weibull  Model) 

The  quantification  of  the  prestressing  effect  can  be  written  as  a 
function. 


i 


I 
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PS  =  f  (O,  'F,  a,  p)  (4) 

where  =  Gauge  Length 
=  Preload  Level 

a  =  Fiber  shape  parameter,  Weibull 
P  =  Fiber  scale  parameter,  Weibull 
The  functional  interrelation  between  these  functions  is  presently 
not  known  and  appropriate  probabilistic  modeling  will  be  required. 

« 
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One  approach  is  to  optimize  each  variable  separately  where  possible. 
For  example,  the  Weibull  strength  parameters  characterize  the  fiber 
strength  distribution.  The  prestress  effect  of  the  preload  level  and  the 
gauge  length  will  alter  these  parameters.  The  optimization  of  the 
effect  on  the  alteration  of  these  parameters  on  the  reliability  of  the 
composite  can  be  investigated  as  one  means  of  optimizing  the  pre- 
stress  effect. 

The  probability  of  failure  (F*)  vs  normalized  load  (x*/P)  graph  for 
the  fibers  (F)  and  composite  (C)  with  the  Failure  Probability  Density 
Function  superimposed  is  shown  in  Figure  5.  This  figure  demonstrates 
how  the  PDF  can  be  transposed  to  a  linear  curve.  The  linear  curve  is 
much  easier  to  use  and  because  it  reveals  in  great  clarity  what  is  hap¬ 
pening  at  the  weak  tail.  This  graph  will  be  used  to  show  the  effect  of 
prestressing  with  respect  to  failure  probability. 


Probability  of  Failure  vs.  Normalized  Strength 
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The  probability  of  failure  for  the  composite  (C)  is  calculated  from 

« 

the  fiber  parameters.  The  composite  reliability  will  diverge  from  the  I 

fiber  reliability  at  Xc  because  of  the  local  load-sharing  stress-concen¬ 
tration  cells.  These  cells  cause  the  composite  fibers  in  the  matrix  to 
break  earlier  because  the  fibers  adjoining  a  broken  fiber  have  to  carry  I 

the  load  carried  by  broken  fiber  by  themselves,  not  spread  equally 
amongst  all  of  the  remaining  fibers.  Until  the  load  Xc,  the  two  curves 
do  not  diverge  and  the  reliability  of  the  composite  is  equal  to  the  reli-  I 

ability  of  the  fibers. 

If  the  fibers  have  been  damaged  so  that  their  distribution  curve  is 
modified  to  appear  like  that  shown  in  Figure  6,  the  result  on  the  com-  j 

posite  strength  can  be  seen.  As  the  composite  graph  swings  upward, 
its  reliability  is  decreasing.  The  composite  is  stronger  than  the  fiber 
due  to  the  load  sharing  but  its  reliability  for  loads  greater  than  Xc  has  j 

diminished.  The  composite  reliability  for  loads  below  Xc  has  not 
changed  at  all. 

If  the  fibers*  strength  distribution  could  be  improved,  the  failure 
probability  curve  would  be  shifted,  as  shown  in  Figure  7.  The  impact 
on  the  composite  strength  is  that  the  curve  swings  downward.  This 
downward  swing  reflects  greatly  increased  reliability.  In  fact,  the 
composite  will  not  fail  theoretically  at  all  for  any  load  less  than  Xc. 

A  method  to  achieve  the  fiber  distribution  in  Figure  7  would  be  to 
proof  test  all  of  the  composite  fibers  up  to  Xc  and  then  remove  any 


25 


Normalized  strength 

Figure  6 


Probability  of  Failure  vs.  Normalized  Strength  (Damaged  Fibers) 


Figure  7 

Probability  of  Failure  vs.  Normalized  Strength  (Improved  Fibers) 
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broken  fibers  from  the  bundle  prior  to  adding  the  matrix.  The  effect  of 
this  process  on  the  failure  probability  curve  is  shown  in  Figure  8  as 
curve  R.  The  cost  of  this  operation,  considering  the  extreme  thinness 
of  the  fibers,  the  susceptibility  to  damage  due  to  handling,  and  the 
extreme  amount  of  fibers  that  make  up  the  composite,  make  this 
process  in  reality  impossible. 

Xp  -  Proof  test  Fibers 
f  _  Original  distribution  f 

R  -  Remove  broken  fibers 
Resultant  distribution 

NRL  -  Not  Remove  broken  fibers 
Long  Gauge  Length 
Resultant  distribution 

NRS  -  Not  Remove  broken  fibers 
Short  Gauge  Length 
Resultant  distribution 

Figure  8 

Effects  on  Failure  Probability  Curves  Due  to  Proof  Testing 

Instead,  what  may  be  done  is  to  preload  the  fibers  to  Xc  on  a  gauge 
length  equal  to  the  ineffective  length  of  the  fiber.  Because  of  the  han¬ 
dling  difficulties,  the  broken  fibers  will  not  be  removed.  The  effect  of 


this  process  on  the  failure  probability  curve  is  shown  in  Figure  8  as 
curve  NRS. 

The  NRS  curve  has  swung  downward  as  compared  to  the  original 
fiber  curve,  but  it  has  not  swung  as  far  down  as  if  the  broken  fibers  had 
been  removed.  The  reason  it  does  not  swing  as  far  is  the  local  load- 
sharing  stress-concentration  cells  that  are  created  when  the  matrix  is 
added. 

If  the  fibers  are  preloaded  to  Xc  but  on  a  gauge  length  that  is 
greater  than  the  ineffective  length,  the  effect  on  the  failure  probability 
curve  is  shown  on  Figure  8  as  curve  NRL.  This  curve  swings  downward 
from  the  original  fiber  curve  even  less  because  of  the  probability  that 
some  fibers  remain  with  segments  that  are  weaker  than  the  preload. 

As  previously  stated,  the  strength  of  the  composite  may  be  calcu¬ 
lated  from  the  strength  of  the  fibers.  If  the  the  fiber  strength  distribu¬ 
tions  are  altered  as  shown  in  Figure  8,  the  composite  strength 
distribution  will  be  altered  in  the  same  direction  as  well.  Using  this,  a 
designer  of  a  composite  may  then  use  a  failure  probability  curve  to 
determine  whether  his  design  meets  his  reliability  specifications.  If  it 
does  not,  he  may  then  decide  to  improve  the  composite  by  prestress¬ 
ing  the  fibers  to  a  critical  load  to  alter  the  failure  probability  to  come 
within  the  requirements. 

The  movement  of  the  failure  probability  curve  can  be  predicted  by 
probability  modeling  based  on  the  preload  level,  the  gauge  length,  the 
ineffective  length,  and  the  original  fiber  strength  parameters.  The 
results  of  such  an  analysis  (by  Wu  and  Harlow)  on  the  fiber  swing  for  a 
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fiber  preloaded  at  the  level  0.3  *  P  .  for  an  a  of  5.0  and  prestressed  at 
multiples  of  its  ineffective  length,  are  shown  in  Figure  9.  The  graph 
demonstrates  that  even  if  the  prestress  gauge  length  is  very  large 
compared  to  the  ineffective  length,  significant  improvement  of  fiber 
strength  distribution  can  be  achieved. 

The  optimization  of  the  prestress  effect  is  to  identify  the  ineffec¬ 
tive  length  or  to  find  the  shortest  prestress  gauge  length  that  can  be 
applied. 
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IV.  EXPERIMENTATION 


The  experimental  work  was  directed  in  two  areas.  The  first  was  to 
measure  the  machine  compliance  of  the  INSTRON  4206  (  a  universal 
materials  testing  machine)  in  order  to  characterize  the  testing 
machine  system  compliance  denoted  by  the  function  8C  (P).  Once  this 
was  completed,  the  second  portion  of  the  experimentation  was  to 
perform  failure  tests  on  AS-4  Graphite  bundle  samples  of  various 
lengths  in  order  to  determine  the  strength  characteristics  and  to 
observe  the  effects  of  gauge  length  on  the  test  results  with  respect  to 
compliance,  slack,  noise,  and  friction.  Bundle  samples  were  tested 
with  no  oil  and  then  with  oil  in  an  attempt  to  reduce  the  effects  of 
friction,  particularly  on  the  longer  samples  of  50  and  25  cm.  Small- 
gauge  lengths  of  2.5  and  1  cm  were  tested  with  particular  interest 
with  respect  to  the  examination  of  their  suitability  to  preloading. 

A.  COMPLIANCE  TESTING 

The  compliance  testing  was  accomplished  using  a  sample  illus¬ 
trated  in  Figure  10.  The  copper  tabs  were  prepared  and  washed  with  a 
dilute  acid  to  enhance  the  performance  of  the  adhesive.  The  graphite 
bundle  was  glued  to  the  lower  set  of  copper  tabs  with  acetylene  adhe¬ 
sive.  The  upper  tabs  were  then  set  so  that  the  distance  (d)  between 
the  the  upper  and  lower  tabs  was  as  near  zero  as  possible.  The 
graphite  bundle  was  then  glued  into  place  with  the  same  adhesive.  The 
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sample  was  prepared  by  Composites  Laboratory  Technician,  Jim 
Nageotte. 


Figure  10 

Compliance  Test  Sample 


The  tensile  test  was  performed  on  the  INSTRON  4206  testing 
|  machine  and  associated  INSTRON  software,  version  4.01  .  A  plot  of  the 

load  vs.  displacement  is  shown  as  Figure  11.  The  data  and  test  sum¬ 
mary  forms  generated  by  the  software  are  enclosed  in  Appendix  E. 

I  The  compliance  data  was  then  fit  to  the  function  shown  as  Equa¬ 

tion  4.  This  was  done  in  accordance  with  the  procedure  for  curve  fit¬ 
ting  detailed  in  Appendix  C.  The  curve  fitting  was  accomplished  with 
I  the  use  of  LOTUS  1-2-3  spreadsheet  program  and  an  IBM/AT.  The 

resulting  coefficients  for  the  curve  fit  of  the  compliance  are  listed  in 
Table  I. 

I 
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Figure  11 

Compliance  Test  Load  vs.  Displacement  Graph 


33 


i 


INSTRON  4206  Compliance  Curve  Fit  Function 


For  0  <  6  <  Dc  1  5  (P)  =  Ai(P)n  +  A2 

For  Del  <  5  5  (P)  =  A3(P)  +  A4  (4) 

TABLE  1 

INSTRON  4206  COMPLIANCE  CURVE  FIT  COEFFICIENTS 


Coefficient 

Value 

Units 

Ai 

.0512 

mm 

A 2 

-.0217 

mm 

n 

1975 

— 

A3 

.0064 

rum /kg 

A4 

.0242 

mm 

B.  BUNDLE  TESTING 

The  bundle  testing  was  accomplished  with  the  use  of  samples  as 
illustrated  in  Figure  12.  The  copper  tabs  were  prepared  and  washed 
with  a  dilute  acid  to  enhance  the  performance  of  the  adhesive.  The 
graphite  bundle  was  glued  to  the  lower  set  of  copper  tabs  with  acety¬ 
lene  adhesive.  The  upper  tabs  were  then  set  so  that  the  distance  (d) 
between  the  the  upper  and  lower  tabs  was  the  target  gauge  length. 
The  graphite  bundle  was  held  straight  with  a  2  kg  weight  while  the 
upper  tabs  were  glued  in  place.  The  samples  were  carefully  stored  in 
the  Composites  Laboratory  until  use.  As  in  the  case  of  the  compliance 
test  sample,  all  of  the  samples  tested  were  prepared  by  Laboratory 
Technician,  Mr.  Jim  Nageotte, 
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The  bundle  tests  were  performed  on  AS-4  Graphite  bundles  of 
3,000  filaments.  The  tests  were  run  in  accordance  with  the  procedure 
outlined  in  Appendix  D.  The  displacement  due  to  compliance  was  sub¬ 
tracted  from  the  total  displacement  and  a  summary  of  the  the  results 
is  listed  in  Table  2.  The  data  and  test  summary  forms  generated  by  the 
software  are  enclosed  in  Appendix  E. 
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TABLE  2 

BUNDLE  TESTING  RESULTS  (COMPLIANCE  REMOVED) 


Material: 

AS-4  Graphite 

(3000  filaments) 

Test  # 

Dry/ 

Gauge 

Modulus 

a 

P 

Comments 

Oil 

Length 

(mm/mm) 

(kg/mm) 

090901 

Diy 

5.0  cm 

1068 

3.214 

.0135 

090902 

Dry 

.5  cm 

* 

090903 

Dry 

.5  cm 

* 

090904 

Dry 

.5  cm 

♦ 

090905 

Dry 

.5  cm 

* 

090906 

Dry 

.5  cm 

* 

090907 

Dry 

5.0  cm 

# 

090908 

Dry 

5.0  cm 

912 

4.465 

.0145 

090910 

Dry 

2.5  cm 

0 

080901 

Oil 

70.0  cm 

Test  Validation 

100901 

Dry 

2.5  cm 

915 

3.000 

.0157 

100902 

Oil 

2.5  cm 

* 

100903 

Oil 

2.5  cm 

* 

100904 

Dry 

50.0  cm 

905 

3.035 

.0077 

100905 

Oil 

50.0  cm 

959 

3.660 

.0085 

100906 

Oil 

50.0  cm 

Not  Analyzed 

100907 

Oil 

25.0  cm 

932 

3.964 

.0101 

100908 

Dry 

25.0  cm 

917 

3.180 

.0097 

*  Adhesive  Failure 


Examination  of  the  load  vs.  percentage  strain  plots  revealed  the 
fact  that  for  the  shorter  gauge  lengths,  2.5  and  I  cm,  the  copper  tabs 
failed  to  grip  the  graphite  bundle  uniformly.  This  caused  some  of  the 
fibers  to  slip  out  of  the  tabs  and  the  data  from  these  tests  were  of  no 
validity.  The  tests  that  had  this  phenomena  are  listed  as  adhesive  fail¬ 
ure  in  the  comments  section  of  Table  2.  Figure  13  shows  a  represen¬ 
tative  example  of  this  phenomenon  as  the  load  remains  constant  while 
the  percent  strain  increases  beyond  the  limit  of  the  experiment. 


Figure  13 

Load  vs.  Percent  Strain  for  Sample  100903  (Adhesive  Failure) 
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An  examination  of  Table  2  shows  that  seven  of  the  18  samples 
resulted  in  satisfactory  tests.  The  seven  successful  tests  results  have 
been  summarized  into  Table  3. 


TABLE  3 

BUNDLE  TESTING  RESULTS  (COMPLIANCE  REMOVED) 

SUCCESSFUL  TESTS 


Material:  AS-4  Graphite  (3000  filaments) 


Test  # 

Dry/ 

Gauge 

Oil 

100901 

Dry 

2.5  cm 

090901 

Dry 

5.0  cm 

090908 

Dry 

5.0  cm 

100908 

Dry 

25.0  cm 

100907 

Oil 

25.0  cm 

100904 

Dry 

50.0  cm 

100905 

Oil 

50.0  cm 

*  Data  not  plotted 


Modulus 

a 

P 

Length 

(kg/mm) 

(mm/mm) 

915 

3.000 

.0157 

1068 

3.214 

.0135 

912 

4.465 

.0145* 

917 

3.180 

.0097 

932 

3.964 

.0101 

905 

3.035 

.0077 

959 

3.660 

.0085 

Figures  14  through  32  are  the  graphs  of  load  vs.  percent  strain  of 
the  test  data,  of  the  compliance,  and  of  the  test  data  with  the  fiber 
removed. 
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AS-4  BUNDLE  TEST  TEST  DATA 


SAkPLE  100901  GAXE  LENGTH  «  15  MU 


Figure  14 

Load  vs.  Percent  Strain  Sample  100901  Test  Data 
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Load  (Kg) 


AS-4  BUNDLE  TEST  COMPLIANCE 


SMuPLE  100901  WGE  LfNGTH  -  25  MM 


Figure  15 

Load  vs.  Percent  Strain  Sample  100901  Compliance 
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AS -4  BUNDLE  TEST  COMPLIANCE  REMOVED 


SAOTHOOM  CAX5C  LENGTH  -  25  UU 


I 


Figure  16 

Load  vs.  Percent  Strain  Sample  100901  Compliance  Removed 
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AS-4  BUNDLE  TEST  TEST  DATA 


SttPlCOHUI  CAJGC  mcm -  so uu 


Figure  17 

Load  vs.  Percent  Strain  Sample  090901  Test  Data 
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Load  (Kq) 


AS-4  BUNDLE  TEST  COMPLIANCE 


SMPLE  090901  G/UGE  LENGTH  -  50  UU 


Figure  18 

Load  vs.  Percent  Strain  Sample  090901  Compliance 
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Load  0<q) 


AS-4  BUNDLE  TEST  COMPLIANCE  REMOVED 


SAMPLE  C9MW  &XJGE  LENGTH  «  50  UU 


Figure  19 

Load  vs.  Percent  Strain  Sample  090901  Compliance  Removed 
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AS— 4  BUNDLE  TEST  TEST  DATA 


SAkPLE  1009 OS  GAUGE  LENGTH  -  150  W 


Figure  20 

Load  vs.  Percent  Strain  Sample  100908  Test  Data 
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AS-4  BUNDLE  TEST  COMPLIANCE 


swpinoogot  w  ibch  .  ik>  i*< 


Figure  21 

Load  vs.  Percent  Strain  Sample  100908  Compliance 
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AS-4  BUNDLE  TEST  COMPUANCE  REMOVED 


$ttPL[10CM»  GAUCC  LOCK  -  250 1*1 


Figure  22 

Load  vs.  Percent  Strain  Sample  100908  Compliance  Removed 
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Load  Ckq> 


AS-4  BUNDLE  TEST  (OIL)  TEST  DATA 


SAUPLC 100907  OK/GC  UK3H  »  20  MU 


Percent  Strain  X 


Figure  23 

Load  vs.  Percent  Strain  Sample  100907  Test  Data 
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Percent  Strain  X 

Figure  24 

Load  vs.  Percent  Strain  Sample  100907  Compliance 
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AS-4  BUNDLE  TEST  (OIL)  COMPLNC  REMOVED 


SAMPLE  10M07  GAUGE  LDCTTH  *  SO  UM 


i  i  l  i  i  l  i  i  i  i  i  i  i  i 


0  0.4  0.6  1.2  1.6  2  2.4 

Percent  Strain  X 

Figure  25 

Load  vs.  Percent  Strain  Sample  100907  Compliance  Removed 
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AS-4  BUNDLE  TEST  TEST  DATA 


SAMPLE  100904  GAUGE  LENGX  -  500  ill 


Figure  26 

Load  vs.  Percent  Strain  Sample  100904  Test  Data 
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Load  (Kg) 


AS-4  BUNDLE  TEST  COMPLIANCE 


SttPLE  100904  GAUGE  IDCH  -  500  kU 


0  0.4  0.i  1.2  l.t  2  2.4 

Percent  Strain  X 


Figure  27 

Load  vs.  Percent  Strain  Sample  100904  Compliance 
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AS -4  BUNDLE  TEST  COMPLIANCE  REMOVED 


SWPLE  100904  GAUGE  LENGTH  -  500  MK 


0  0.4  Q.i  U  U  1  2.4 


Percent  Strain  X 


Figure  28 

Load  vs.  Percent  Strain  Sample  100904  Compliance  Removed 


<D>0  pooh 


AS-4  BUNDLE  TEST  (OIL)  TEST  DATA 


SAAPLE  100905  GAUGE  LOCK  -  500  IN 


Figure  29 

Load  vs.  Percent  Strain  Sample  100905  Test  Data 
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AS— 4  BUNDLE  TEST  (OIL)  COMPLIANCE 

SAiPtI  100305  GAUGE  LOCH  -  500  Mi 

21- 

24- 

22- 

20- 

15  - 

IB  - 
a 

TJ 

o  u- 

j 

10- 
5- 
B- 
4- 
2- 
0- 

0  0.4  0.5  1.2  1.B  2  2.4 

Percent  Srcin  % 


/ 
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Figure  30 

Load  vs.  Percent  Strain  Sample  100905  Compliance 
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AS-4  BUNDLE  TEST  (OIL)  COMPLNC  REMOVED 


SAMPLE  10905  GAUGE  LENGTH  ■  5D0  MM 


0  0.4  0.5  1.2  1.5  2  2.4 

Pferctnt  Strain  X 


Figure  3 1 

Load  vs.  Percent  Strain  Sample  100905  Compliance  Removed 
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V.  CONCLUSIONS 
A.  COMPLIANCE  TESTING 

The  curve-fitting  procedure  outlined  in  Appendix  C  was  successful 
in  arriving  at  a  continuous  function  to  describe  the  load-induced  dis¬ 
placement  of  the  load  train  for  the  INSTRON  testing  machine.  Figure 
32  shows  the  data  points  measured  during  the  zero  gauge  length  test. 
Figure  33  shows  the  curve  fit  function  using  the  parameters  listed  in 
Table  2.  Figure  34  then  shows  the  data  points  and  the  curve  fit 
superimposed  on  one  another.  The  greatest  deviation  in  load  is 
.0017  kg. 

Based  on  the  small  deviation  of  the  data  points  to  the  curve  fit 
function,  the  curve  fit  of  the  compliance  is  satisfactory  to  be  used  to 
subtract  the  displacement  due  to  compliance  from  bundle  test  data. 

The  curve  fit  solution  for  the  machine  compliance  is  satisfactory 
for  loads  greater  than  .1  kg.  The  curve  fit  does  not  do  a  good  job  in  the 
region  below  this  load.  This  is  due  to  the  inconsistent  zero  displace¬ 
ment  length.  The  zero  point  for  displacement  is  very  difficult  to 
achieve  for  bundle  testing.  The  sensitivity  of  the  load  cell  is  insuffi¬ 
cient  to  consistently  identify  an  actual  zero  load  level  on  the  bundle. 
The  effect  of  this  on  the  fiber  strength  distribution  parameters  is  that 
some  elements  of  the  lower  tail  are  distorted. 
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Figure  33 

INSTRON  Compliance  Curve  Fit  Function 
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Figure  34 

0  Gauge  Length  INSTRON  Compliance  Test  and  Curve  Fit 
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The  compliance  testing  method  for  the  INSTRON  has  been 
demonstrated  to  be  feasible  and  to  provide  realistic  results.  Testing 
should  continue  to  determine  whether  the  compliance  is  repeatable  or 
whether  it  is  a  function  of  other  factors  such  as  temperature,  humid¬ 
ity,  or  age  of  the  machine.  Full  understanding  of  compliance  is  impor¬ 
tant  because  of  the  magnitude  of  compliance  effects  on  short  gauge 
lengths;  it  is  the  short  gauge  lengths  to  which  stressing  is  applied. 

B.  BUNDLE  TESTING 

The  bundle  tests  demonstrated  the  difficulty  in  obtaining  test  data 
for  short  gauge  lengths.  Table  2  illustrated  the  fact  that  no  successful 
1  cm  tests  were  completed  and  only  one  2.5  cm  test  was  completed. 
An  improved  technique  of  sample  preparation  to  prevent  the  slipping 
of  the  fibers  through  the  tabs  is  needed.  Possibly  an  adhesive  with  bet¬ 
ter  wetting  qualities  could  be  tried. 

The  difficulty  in  testing  the  small  gauge  lengths  highlights  two 
important  effects  on  prestressing.  The  first  is  that  the  shorter  the 
prestress  gauge  length,  the  more  difficult  the  process  becomes  due  to 
the  difficulty  in  maintaining  gauge  length.  The  second  important  effect 
is  that  the  prestress  effect  on  fiber  parameter  improvement  will  have 
to  be  based  on  a  gauge  length  less  than  the  ineffective  length.  The 
optimization  process  will  now  Include  another  variable,  the  minimum 
possible  gauge  length  due  to  machine  constraints.  As  the  minimum 
gauge  length,  increases  the  effectiveness  of  the  prestress  will 
decrease. 
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The  Importance  of  accurately  removing  the  displacement  in  the  I 

test  data  for  short  gauge  lengths  was  demonstrated  by  the  bundle  •( 

tests.  Figures  15.  18.  21,  24.  27,  and  30  show  that  as  a  function  of  1 

percent  strain,  the  shorter  the  gauge  length,  the  greater  the  percent-  I 

age  of  test  displacement  is  due  to  machine  compliance.  Figure  15  ■ 

shows  that  the  percent  strain  displacement  due  to  compliance  for  the 
shortest  gauge  length  (2.5  cm)  was  greater  than  .6  percent,  or  nearly 
25  percent  of  the  total  test  displacement.  Figure  27  shows  that  the  JJ 

percent  strain  displacement  due  to  compliance  for  the  longest  gauge 
length  (50  cm)  was  less  than  .02  percent,  which  is  less  than  1  percent 
of  the  total  test  displacement.  fl 

The  inaccuracy  incurred  in  not  eliminating  the  compliance  dis¬ 
placement  for  long  samples  before  reducing  the  data  to  determine  the 
fiber  strength  parameters  is  minimal.  In  the  case  of  short  samples,  • 

however,  this  is  not  the  case.  If  the  compliance  displacement  is  not 
removed,  the  effect  on  the  strength  parameters  will  be  the  same  as  if 
the  gauge  length  was  actually  longer  than  what  was  tested.  This  results  1 

in  the  fiber  parameters  being  conservative  in  the  prediction  of  fiber 
strength.  While  this  is  safe,  it  defeats  the  purpose  of  prestressing. 

Figures  25  and  31  are  the  graphs  of  the  oiled  bundles.  Comparing  • 

these  graphs  to  Figures  22  and  28,  which  are  the  dry  or  not  oiled 
samples,  the  effects  of  friction  can  be  seen  in  the  shape  of  the  curves 
after  the  maximum  load  has  been  reached.  The  dry  sample  load  drops  jj 

in  large  vertical  jumps.  These  large  jumps  correspond  to  many  fibers 
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breaking  at  once  instead  of  sequentially.  Figure  35  shows  the 
simulated  test  data  for  a  100  mm  sample.  The  curve  is  very  smooth 
and  does  not  show  the  sharp  vertical  drops  the  actual  test  data  shows. 

The  large  vertical  drops  in  the  actual  test  data  are  noticeably 
reduced,  but  not  eliminated,  in  the  oiled  samples.  The  oil,  a  soybean 
derivative,  reduces  the  friction  but  does  not  eliminate  all  of  the  effects. 
The  effect  of  friction  on  the  test  does  appear  to  dominate  the  tail  end 
or  the  final  30  percent  of  the  percent  strain.  The  effect  of  friction  on 
the  fiber  strength  parameters  is  difficult  to  predict.  One  effect  that  is 
clear  is  that  the  data  for  the  first  50  percent  to  60  percent  of  the 
breaks  in  the  long  samples  do  not  appear  to  be  affected  by  friction  as 
the  remaining  breaks.  In  addition,  Figures  16  and  19  indicate  that  the 
effect  of  friction  declines  as  the  gauge  length  becomes  shorter. 

The  need  to  eliminate  friction  between  the  fibers  in  a  bundle  test 
is  important  if  long-gauge  length  tests  are  to  be  accurately  completed. 
Different  oils  may  be  tried  as  lubricants  to  improve  upon  the  perfor¬ 
mance  of  the  soybean  oil.  If  the  friction  cannot  be  eliminated,  it  must 
be  brought  down  to  a  minimal  level. 

The  effect  of  noise  in  the  bundle  tests  appears  to  be  minimal.  Fig¬ 
ures  36,  37,  and  38  are  graphs  of  simulated  bundle  tests  with  noise 
added  to  the  load  data.  Figure  4  is  a  graph  of  the  same  simulation  with 
no  noise  present. 

A  noise  level  of  1  gm  is  equivalent  to  .002  percent  for  the  50  kg 
load  cell  used  by  the  INSTRON.  The  difference  in  the  curves  from 
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Figure  36 

Bundle  Test  Simulation  With  ±1  gm  Noise 
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Figure  37 

Bundle  Test  Simulation  With  ±10  gm  Noise 
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BUNDLE  SIMULATION  DISCRETE  DATA 


Figure  38 

Bundle  Test  Simulation  With  ±100  gm  Noise 
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DISPLACEMENT  (M^) 


Figure  35  to  Figure  36  is  barely  noticeable.  A  noise  level  of  10  gms  is 
equivalent  to  .02  percent.  Figure  37  shows  some  effects  of  noise,  but 
they  are  minimal.  A  noise  level  of  100  gms  is  equivalent  to  .2  percent. 
This  is  an  excessive  amount  of  noise  and  the  effects  of  a  noise  level 
this  high  are  clearly  shown  by  Figure  8. 

The  INSTRON  manual  specifies  the  sensitivity,  or  noise  level,  of 
the  load  cell  to  be  .01  kg.  This  corresponds  to  roughly  the  noise  level 
of  Figure  37.  An  inspection  of  Figure  37  demonstrates  that  the  effect 
of  this  noise  level  is  barely  discemable  when  compared  to  Figure  35,  a 
simulation  with  no  noise. 

This  research  has  demonstrated  that  the  physical  process  of 
preloading  fibers  with  a  short  gauge  length  is  difficult  due  to  the 
physical  limitations  of  the  sample  gripping  mechanism.  More  work  is 
required  to  find  an  adequate  means  to  grip  or  restrain  the  fibers  for 
short  gauge  lengths. 

The  removal  of  the  displacement  due  to  compliance  in  the  test 
results  is  significant  for  gauge  lengths  shorter  than  2.5  cm  in  order  to 
determine  accurately  the  fiber  strength  distribution  parameters.  The 
effect  of  compliance  diminishes  with  longer  gauge  lengths.  The  fric¬ 
tion  between  the  fibers  in  bundle  alters  the  shape  of  the  load  vs.  dis¬ 
placement  curve.  Friction  has  more  of  an  effect  on  longer  samples 
than  on  short  samples.  Friction  may  be  reduced  by  using  oil  on  the 
samples,  but  it  does  not  entirely  eliminate  the  effect.  Noise  is  not  a 
significant  factor  in  altering  the  data  and  does  not  interfere  with  the 
strength  distribution  parameters. 
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VI.  RECOMMENDATIONS 

Further  study  into  the  magnitude  of  the  machine  compliance  for 
the  INSTRON  4206  testing  machine  should  be  initiated.  The  effects  of 
temperature,  the  gripping  devices,  and  repeatability  should  be  investi¬ 
gated.  To  complete  this,  a  quicker  method  to  determine  the  curve  fit 
coefficients  needs  to  be  developed.  A  Fortran  program  that  can  accu¬ 
rately  and  quickly  solve  the  curve -fitting  routine  in  Appendix  C  should 
be  written  and  validated.  A  series  of  compliance  tests  may  then  be 
completed  to  build  a  database. 

The  methods  for  sample  preparation  and  bundle  testing  need  to 
become  less  technician  dependent.  There  are  too  many  areas  where 
mishandling  of  the  samples  may  introduce  errors  into  the  data  by 
damaging  the  fibers.  The  reliance  on  the  laboratory  technician  to  build 
the  samples,  align  the  samples,  and  then  set  the  gauge  length  all  add 
variability  to  the  data  that  should  be  reduced  if  not  eliminated. 
Hydraulic  grips  can  eliminate  the  hand-  and  wrench-tightening  pro¬ 
cedure  required  by  the  current  grips.  A  brace  or  stand  should  be 
designed  and  built  to  hold  the  sample  during  loading  into  the 
INSTRON.  This  will  eliminate  much  of  the  handling  the  sample  must 
now  endure. 

The  prestress  effects  of  preloading  the  fibers  to  a  load  greater 
than  the  critical  load  should  be  investigated.  This  is  a  probability 
problem  whose  solution  could  be  integrated  into  the  bundle  simulation 
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program.  This  program  could  then  be  used  to  predict  the  prestressing 
performance  and  then  validated  with  tests  performed  in  the  laboratory 
on  gauge  lengths  that  can  be  successfully  tested.  This  is  very  important 
because  the  poor  test  results  of  the  short  samples  indicate  that  a  short 
preload  gauge  length  (<  1  cm)  may  be  difficult  to  achieve. 

The  effects  of  prestressing  with  respect  to  long-term  reliability 
and  composite  fatigue  failure  should  be  investigated.  The  impact  of 
local  load  stress  cells  on  long-term  reliability  needs  to  be  addressed 
quickly.  The  urgency  of  this  aspect  of  prestressing  is  that  due  to  the 
nature  of  fatigue  failure,  data  taking  is  a  very  long  process. 
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APPENDIX  A 


MODEL  SELECTION 

Composite  fiber  materials  such  as  graphite,  aramid,  or  boron  are 
characteristically  long  and  thin  when  compared  to  the  dimensions  of 
the  composite.  For  example,  the  diameter  of  a  graphite  fiber  is  on  the 
order  of  5  microns,  while  the  fiber  length  can  extend  to  hundreds  of 
feet  for  a  1  square  foot  pressure  vessel.  Under  tension,  the  fibers  fail, 
or  break,  at  positions  where  flaws  or  imperfections  exist.  These  flaws 
cause  the  fiber  to  break  at  loads  much  below  the  theoretical  strength 
of  the  fiber,  which  is  the  chemical-bonding  strength.  The  type  of  the 
flaw  may  be  scratching,  chemical  erosion,  contamination,  or  one  of 
many  others.  The  source  of  the  flaw  may  have  come  in  production, 
handling,  or  service  use.  The  spatial  location,  density,  and  severity  of 
these  flaws  determine  the  strength  of  each  fiber. 

A  single  fiber  or  filament  may  be  visualized  as  a  chain  made  up  of 
many  links.  These  links,  or  segments,  each  have  a  unique  strength 
that  is  limited  by  the  severity  of  the  flaws  contained  in  that  segment. 
The  more  severe  the  flaw,  the  weaker  the  segment.  For  a  long  fiber, 
one  with  many  segments,  the  strength  of  that  fiber  is  only  the 
strength  of  the  most  severe  flaw,  or  the  weakest  link.  If  the  severity  of 
the  flaw  is  known  to  fall  within  a  certain  probability  distribution  and 
the  location  of  the  flaws  are  randomly  distributed  along  the  fiber,  then 
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a  “weak  link"  model  can  be  assumed  for  predicting  the  strength  of  a 
fiber. 

From  statistics,  the  Weibull  distribution  best  fits  the  “weak  link" 
model.  Phoenix  and  Wu,  in  Reference  2,  describe  the  Weibull  statisti¬ 
cal  model  as  it  applies  to  composite  fiber  strength  prediction  in  detail. 
In  general,  the  fiber  is  partitioned  into  a  series  of  segments  (m).  The 
shortest  limiting  length  segment  is  equal  to  the  ineffective  length  or 
effective  load  transfer  length  of  the  fiber.  This  is  on  the  order  of  mm 
for  graphite/epoxy.  A  bundle  is  then  formed  by  joining  several,  (n), 
fibers  in  parallel.  For  example,  in  this  project.  AS-4  bundles  where  m 
equals  i  and  n  equals  3,000  were  tested  in  the  laboratory. 

Based  on  the  Weibull  weak  link  model,  the  failure  cumulative  den¬ 
sity  function  (CDF)  of  a  fiber  can  be  explicitly  written  as  a  function  of 
the  two  Weibull  parameters  a  and  ft.  The  Weibull  failure  CDF  is  written 
as  Equation  5. 


F  =  1  -  exp  {-(x/p)a) 


(5) 


The  two  parameters,  a  and  p,  correspond  to  the  distribution 
shape  and  scale,  a  adjusts  the  shape  of  the  distribution  by  altering  the 
skew  and  range  of  values.  As  shown  in  Figure  39,  a  high  shape  para¬ 
meter,  or  a,  corresponds  to  a  very  narrow  distribution;  an  a  equal  to 
3.5  approximates  a  normal  distribution;  and  a  smaller  a  corresponds 
to  a  veiy  wide  distribution  with  large  tails.  When  a  is  large,  the  scale 
parameter  (3  is  approximately  the  mean.  The  AS-4  that  was  tested  in 
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Figure  39 

Effect  of  the  Shape  Parameter  on  the  Weibull  Distribution 

The  Weibull  distribution  is  not  only  a  physically  appropriate  model 
for  fiber  strength,  it  has  an  additional  advantage  in  that  its  failure 
cumulative  distribution  (CDF)  can  be  written  as  an  explicit  equation.  In 
addition,  the  CDF  may  be  linearized  with  respect  to  p.  This  results  in 
Equations  6  and  7.  These  two  equations  may  then  be  graphed  together 
to  form  a  graph  such  as  the  one  shown  as  Figure  40.  This  graph  of  the 
F*/x*  function  allows  direct  graphical  correlation  of  strength  to  failure 
probability. 


Figure  40 

Weibull  Linearized  Failure  Prediction  Graph 


F*  =  ln(-ln(l  -  exp  {-(P/p)“))) 

(6) 

x*  =  In  (P/p) 

(7) 

where  P  is  the  load  subjected  to  the  filament. 

To  predict  reliability,  which  is  one  minus  the  failure  probability 
for  a  fiber,  the  a  and  J3  must  be  known  for  that  fiber  to  a  sufficient 
degree  of  accuracy.  Fiber  bundle  failure  tests  are  performed  in  the  lab¬ 
oratory  to  arrive  at  the  parameters.  A  difficulty  arrives  with  the  fact 
that  the  parameters  are  a  function  of  the  length  of  a  fiber.  This  corre¬ 
sponds  to  the  fact  that  a  long  chain  is  weaker  than  a  short  chain 
because  the  longer  chain  has  a  higher  probability  of  having  a  weak 
segment. 
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The  Weibull  model  provides  a  means  to  normalize  the  parameters 
between  different  lengths  of  a  common  fiber  if  the  parameters  are 
known  for  any  one  length.  This  allows  for  the  laboratory  results  to  be 
applied  to  fibers  of  lengths  different  than  the  lengths  tested.  This  fea¬ 
ture  of  the  Weibull  function  is  commonly  called  the  size  effect. 

If  the  strength  distribution  parameters  (11  and  ocl  are  known  for  a 
fiber  of  length  li  and  the  parameters  (I2  and  c*2  are  desired  for  the 
same  fiber  but  of  a  different  length  I2  the  following  derivation  provides 
the  relationship: 


Fi  =  1  -  exp  {-(x/pi)ai) 

(5) 

R  =  1  -  F 

(Reliability  =  1  -  Failure) 

(8) 

Ri  =  exp  {-(x/pi)ai) 

(9) 

m  =  I1/I2 

(10) 

Rt  =  flRi  R2 . =  (Ri)m 

Rt  =  (exp  {-(x/Pi)ai)]m 

(ID 

R2  =  exp  {-m(x/pi)ai} 

(12) 

R2  =  exp  {-(x/p2)a2) 

(9) 

Equating  Equations  9  and  12  then: 
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exp  {-m(x/Pi)ai}  =  exp  {-(x/^2>a2} 
{-m(x/pi)«i}  =  {'(x/p2)a2} 
(x/Cm)-1^!  p!)«i  =  (x/p2  )a2 


(13) 

Solving  Equation  13  results  in  Equations  14  and  15. 


(X2  =  ai 

(14) 

p2  =  Pi(m)-1/ai 

(15) 

An  implication  of  the  size  effect  relations,  Equations  14  and  15,  is 
that  a  long  composite  made  up  of  many  segments  may  be  reduced  to  a 
composite  of  one  segment.  This  simplification  provides  for  easier 
mathematical  modeling  of  the  composite  as  a  large  composite  model 
shown  in  Figure  41  with  many  segments  may  be  reduced  to  the  much 
smaller  model  shown  in  Figure  42.  This  greatly  reduces  the  amount  of 
computational  effort  required  to  model  composite  strengtn  or  bundle 
strength  as  the  model  matrix  has  been  reduced  from  a  (m  X  n)  to  a  (1 
X  n). 
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Figure  42 

Composite  Model  With  Size  Effect 
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APPENDIX  B 


BUNDLE  FAILURE  SIMULATION 

The  purpose  of  the  bundle  failure  simulation  was  twofold.  The  first 
was  to  allow  the  user  to  readily  produce  simulated  bundle  test  data 
with  variable  fiber  strength  distribution  parameters.  The  resulting  data 
could  then  be  graphically  presented  into  load  versus  displacement  or 
load  verses  percent  strain  representations.  These  graphs  were  utilized 
to  gain  understanding  and  insight  into  the  effect  of  fiber  strength 
parameters  and  gauge  length  on  fiber  failure  process. 

The  second  purpose  of  the  bundle  failure  simulation  was  to 
provide  a  means  of  verification  for  data  interpretation  and  reduction 
software  being  developed  for  use  in  conjunction  with  the  INSTRON 
Testing  machine  used  in  the  Composites  Laboratory,  Naval 
Postgraduate  School.  The  simulation  could  be  used  to  ensure  that  the 
software  returned  the  parameters  from  the  simulated  test  data  that 
were  used  as  inputs  to  the  bundle  simulation.  For  this  project,  the 
compliance  removal  procedures  and  software  described  in  Appendices 
C  and  D  were  validated  with  this  program.  Lt.  Joseph  Schmidt,  in  a 
separate  project,  used  this  program  to  validate  software  that  reduced 
bundle  test  data  into  the  fiber  strength  probability  distribution 
parameters. 

The  bundle  failure  simulation  program  “T2"  takes  fiber  strength 
distribution  parameters  a  and  p  for  the  length  of  the  bundle,  the 
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number  of  filaments,  the  length  of  the  bundle,  the  fiber  modulus,  and 
testing  machine  settings  of  cross-head  speed  and  sampling  rate,  and 
then  produces  two  sets  of  load  vs.  displacement  data.  The  first  set  is 
the  loaded  and  displacement  determined  for  each  individual  filament 
break.  This  is  referred  to  as  the  analog  data.  The  second  set,  which 
simulates  fiber  testing  machine  data,  is  the  continuous  load  and 
incremental  displacement  data,  which  is  referred  to  as  discrete  data. 
The  program  also  has  the  added  feature  of  inputting  a  noise  range.  The 
maximum  number  of  filaments  is  1,000. 

The  program  solves  for  the  strength  of  each  filament  based  on  the 
Weibull  distribution  and  the  size  effect.  An  explanation  of  the  Weibull 
function  and  of  the  size  effect  is  detailed  in  Appendix  A. 

A.  ANALOG  DATA  SIMULATION 

The  analog  data  was  determined  by  using  a  random  number 
generator  to  to  generate  a  random  strength  of  each  fiber,  determine 
the  displacement  of  each  break,  and  then  compute  the  load  of  the 
bundle  at  each  break. 

Equation  16  was  used  to  determine  the  strength  of  each  fiber. 

Pf  =  ln(-ln(l  -  exp  {-(x/P)a)))  (16) 

where  Pf  =  failure  load  of  the  the  fiber 
x  =  random  number  0  <  x  <  1 

Equation  17,  which  is  a  manipulation  of  Hooke’s  law,  is  used  to 
determine  the  displacement  of  each  fiber  at  the  failure  load. 
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df  =  Pf  *  1/E 


(17) 


where  df  =  Displacement  of  the  Fiber 

1  =  Length  of  the  Fiber  (Gauge  Length) 

E  =  Strain  Modulus 

The  displacement  at  each  break  was  stored  in  an  array  and  then 
sorted  in  ascending  order.  The  load  of  the  bundle  at  each  break  was 
then  determined  by  using  Equation  18.  This  is  also  a  manipulation  of 
Hooke's  Law.  At  each  break,  because  the  simulation  is  for  a  controlled 
displacement,  the  load  in  the  bundle  drops  by  a  factor  of  1/n.  This 
data  point  is  solved  for  by  Equation  19. 

(Pbh  =  ((df)i  *  E  *  (n-i  +1))/1  (18) 

(Pbh+  =  ((dfh*  E  *  (n-i))/l  (19) 

where  i  =  1,2,3,  .  .  .  n  (Number  of  Breaks) 

B.  DISCRETE  DATA  SIMULATION 

The  discrete  data  was  found  by  utilizing  the  analog  data.  The 
analog  data  was  transformed  into  a  continuous  function  by  connecting 
the  data  points.  The  displacement  between  data  points  is  found  by 
Equation  20. 

Deld  =  CHS/SR  (20) 
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where  CHS  =  Cross-head  Speed  (mm /min) 

SR  =  Sampling  Rate  (pts/sec) 

The  discrete  data  displacement  is  in  increments  of  Deld  while  the 
discrete  data  load  is  the  intercept  of  the  analog  continuous  data  for 
that  appropriate  displacement.  Equation  21  is  used  to  compute  the 
discrete  displacement  points.  Discrete  data  points  are  solved  for  until 
the  load  reaches  zero. 

d  =  Deld  *  I  1  =  0.1.2 _ Load  =  0  (21) 

C.  NOISE  SIMULATION 

The  user  has  the  option  of  simulating  the  effects  of  noise  in  the 
load  transducer  of  a  bundle  test.  If  desired,  the  noise  is  added  to  the 
discrete  load  data  by  Equation  22. 

P noise  =  P  -  ((Ph  '  Pi)/ 2)  +  X  *  (Ph  *  Pi)  (22) 

where  P  =  Original  Load 

Ph  =  High  Noise  Range 
Pj  =  Low  Noise  Range 

D.  PROGRAM  LISTING 

The  program  was  written  in  Microsoft  Fortran  4.01.  It  is  designed 
to  be  operated  from  any  IBM/Compatible  computer  with  two  floppy 
drives  or  hard  disk  and  one  floppy  drive.  The  executable  program  is 
run  from  the  A:  drive  or  the  hard  disk  and  the  data  is  written  to  the  B: 
drive.  To  execute  the  program  type:  T5  <enter>. 
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SUBROUTINE  LISTING: 


INIT 

Initializes  arrays  to  zero 

INPUT 

Prompts  user  to  input  required  value 

RAND 

Random  number  generator  (user  must  input  seed) 

SLOAD 

Solves  for  the  breaking  stress  of  each  filament 

SORT 

Sorts  the  filaments  by  breaking  stress  and  then 
solves  for  breaking  displacement 

ANLDAT 

Creates  the  analog  load  and  displacement  file 

DISDAT 

Solves  for  the  continuous  and  creates  the  discrete 
load  and  displacement  file 

NOISE 

Allows  user  to  input  random  noise  to  a  specified 
threshold  into  the  discrete  load  data 

OUTPUT 

Writes  the  analog  and  discrete  data  to  floppy 
disks 

VARIABLE  LISTING: 

DATA(0: 1000,7) 

Workspace  array 

ADATA(0:2000,2) 

Analog  data  output  file;  Displacement,  Load,  listed 
sequentially  by  first  break  to  last 

DDATA(0:2000,2) 

Discrete  data  output  file;  Displacement,  Load, 
listed  as  a  function  of  sampling  time 

ALPHA 

Weibull  shape  parameter  for  the  given  length 

BETA 

Weibull  scale  parameter  for  the  given  length 

E 

Modulus  of  fiber 

L 

Length  of  the  Fiber 

N 

CHS 

Number  of  filaments  in  the  bundle  (max,  =  1,000) 

Simulated  testing  machine  cross-head  speed 

CHS 


SR 

DELD 

POINTS 

SEED 

SEED1 

SEED2 

BREAKS 

LIMIT 

PH 

PL 

OUTPUT  FILES 

ADATA 

DDATA 


Simulated  machine  sampling  rate 

Distance  between  data  points  for  discrete  data 

Number  of  discrete  data  points 

Random  number  generator  seed  for  strength 
computations 

Random  number  generator  seed  for  discrete  data 
load  during  break  generator 

Random  number  generator  seed  for  noise 
generator 

Counter  that  is  used  to  find  slope  of  analog  data 
for  up  to  n- 1  breaks 

Maximum  strength  of  fiber  i  used  in  solving  for 
the  maximum  displacement  for  n-1  breaks 

Upper  noise  range 

Lower  noise  range 


Analog  data;  displacement,  load 
Discrete  data;  index,  displacement,  load 
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PROGRAM  T2 


*  PROGRAM  T2  WILL  SIMULATE  A  BUNDLE  TEST  FOR  N  FIBERS  IN  PARALEL  * 

*  FOR  GIVEN  FIBER  STRENGTH  PARAMETERS  AND  LENGTH  OF  THE  BUNDLE. 

*  THE  PROGRAM  WILL  OUTPUT  ANALOG  DATA  THAT  IS  A  LISTING  OF  THE 

*  DISPLACEMENT  AND  LOAD  OF  EACH  FIBER  BREAK.  THE  PROGRAM  WILL 

*  OUTPUT  DISCRETE  DATA  THAT  SIMULATES  MACHINE  TESTING  DATA  FOR  THE  * 


*  BUNDLE.  THE  OUTPUT  IS  BY  INDEX,  DISPLACEMENT  AND  LOAD.  THE 

*  NUMBER  OF  DATA  POINTS  IS  A  FUNCTION  OF  THE  SIMULATED  CR03SIIEAD 

*  SPEED  AND  THE  SAMPLING  RATE.  * 

*  * 

*  THE  OUTPUT  IS  WRITTEN  TO  DRIVE  B:  ADATA.  AND  DDATA.  * 

*  * 


***************************** 


*  ******  **  *  *  *  **  *  V?  ***********  ******* 


*  DECLARATION  STATEMENTS 

REAL  DATA(0; 1000,7) ,  DDATA(0;  5000,2) .ALPHA, BETA, L, CHS, PH, PL, 
IE ,DELD, SR,ADATA( 0; 2000,2) 

INTEGER  N, I, J,K, POINTS, Q 
1  CONTINUE 


*  PRINT  BANNER  AND  START  PROGRAM 

PRINT-'-, ’COMPOSITE  BUNDLE  TEST  SIMULATION  PROGRAM' 

CALL  INIT(DATA, DDATA, ADATA) 

CALL  INPUT! ALPHA , BETA , L , CHS . E , N , SR) 

CALL  RAND ( DATA, N) 

CALL  S LOAD (DATA , N , ALPHA , BETA) 

CALL  SORT( DATA, N,E,L) 

CALL  ANLDAT! DATA ,N) 

CALL  DI SDAT! DDATA , DATA , N , E , L , CHS , SR , POINTS ) 

CALL  NOISE! DDATA, POINTS) 

CALL  OUTPUT! DATA , DDATA ,N,E,L, ALPHA , BETA , POINTS , ADATA ) 

*  ALLOW  USER  TO  RERUN  PROGRAM 

PRINT* 

PRINT*,’ RUN  COMPLETE,  DO  YOU  WISH  TO  RUN  ANOTHER  SET  OF  DATA’ 
PRINT*,’  IF  YES  TYPE  1,  ELSE  TYPE  0  TO  EXIT  PROGRAM’ 

READ*,  Q 

IF( Q  . EQ. 1)  GOTO  1 

*  PRINT  BANNER  AND  END  PROGRAM 

PRINT* 

PRINT*,'  END  OF  RUN,  HAVE  A  NICE  DAY' 

END 

*  SUBROUTINE  INIT 

SUBROUTINE  INIT  INITIALIZES  THE  ARRAYS  TO  ZERO 


*  DECLARATION  STATEMENTS 

SUBROUTINE  INIT! DATA, DDATA, ADATA) 

REAL  DATA! 0; 1000 , 7 ), DDATA! 0; 5000 , 2) , ADATA! 0; 2000,2) 
INTEGER  I.J 


PRINT  BANNER  AND  START  SUBROUTINE 


PRINT* * * * * 

PRINT*,’*****  RUNNING  INITIALIZATION  SUBROUTINE  ****** 

PRINT* 

DO  10  I  =  0,1000 
DO  9  J  =  1,7 
DATA( I , J)  =  0. 

9  CONTINUE 
10  CONTINUE 

DO  20  I  =  0,  5000 

DO  19  J  =  1,2 

DDATAC I , J)  =  0. 

19  CONTINUE 

20  CONTINUE 

DO  30  I  =  0,  2000 

DO  29  J  =  1,2 

ADATA( I , J)  =  0 

29  CONTINUE 

30  CONTINUE 
END 
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SUBROUTINE  INPUT 

*  SUBROUTINE  INPUT  PROMPTS  THE  USER  TO  INPUT  THE  REQUIRED 

*  DATA  TO  PERFORM  THE  SIMULATION.  INPUTS  ARE: 


ALPHA,  BETA,  N,  L,  E  -  FIBER  PARAMETERS 

CHS,  SR  -  TESTING  MACHINE  PARAMETERS 


*  DECLARATION  STATEMENTS 

SUBROUTINE  INPUT! ALPHA , BETA , L , CHS , E , N , SR ) 

REAL  ALPHA, BETA, CHS, L, SR, E 
INTEGER  N 

*  PRINT  BANNER  AND  START  SUBROUTINE 

PRINT* 

PRINT*, ******  RUNNING  INPUT  SUBROUTINE  ****** 

PRINT* 

PRINT*,' INPUT  ALPHA  AND  BETA  FIBER  PARAMETERS' 

READ*,  ALPHA,  BETA 

PRINT*, ' INPUT  NUMBER  OF  FILAMENTS  IN  BUNDLE  (MAX  =  1000)’ 
READ*,  N 

PRINT*, ’ INPUT  MEAN  BUNDLE  LENGTH  (CM)’ 

READ*,  L 
L  =  L*10 

FRINT*, ' INPUT  LOAD  MODULUS  OF  FIBER  (GM/MM/MM)’ 

READ*,  E 

PRINT*, 'INPUT  SIMULATION  CROSSHEAD  SPEED  (MM/MIN) ' 

READ*,  CHS 

PRINT*, ' INPUT  SIMULATION  SAMPLING  RATE  (PTS/SEC)' 

READ*,  SR 
PRINT* 

PRINT* , *  INPUT  COMPLETE  * 

PRINT* 


I 


END 
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*  SUBROUTINE  RAND 

* 


* 

* 


*  SUBROUTINE  RAND  GENERATES  A  RANDOM  NUMBER  FROM  0  TO  1  FOR  EACH  * 

*  FILAMENT  IN  THE  BUNDLE.  THE  RANDOM  NUMBER  GENERATOR  IS  STARTED  * 

*  FROM  THE  INITIAL  INPUT  SEED.  THE  FUNCTION  RAND  IS  USED  TO  * 

*  ACTUALLY  FIND  THE  RANDOM  NUMBER.  * 

*  * 


**********************Vr.Y*****'ftyr**iY********yr*>YiV****A*iY,>ViV*yr*************Vc 

SUBROUTINE  RAND ( DATA, N) 


*  DECLARATION  STATEMENTS 

REAL  DATA(0; 1000, 7) ,RNG 
INTEGER  N,  SEED,  I 


*  PRINT  BANNER  AND  START  SUBROUTINE 
PRINT* 

PRINT*,'*****  RUNNING  RANDOM  NUMBER  GENERATOR  SUBROUTINE  ****** 
PRINT* 

PRINT* 

PRINT* , ' INPUT  ANY  ODD  INTEGER  (STRENGTH  RNG  SEED)' 

READ*,  SEED 
PRINT* 

DO  100  I  =  1,  N 
DATA(I.l)  =  I 
DATA(I,3)  =  RNG( SEED) 

100  CONTINUE 


END 
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*  SUBROUTINE  SLOAD 
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Vr 


*  SUBROUTINE  SLOAD  WILL  GENERATE  A  BREAKING  STRENGTH  FOR  EACH 

*  FILAMENT  AND  FROM  THE  MODULUS  FIND  THE  CORESPONDING  DISPLACEMrNT.  * 
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SUBROUTINE  SLOAD (DATA,  N, ALPHA, BETA) 


*  DECLARATION  STATEMENTS 

REAL  DATA(0; 1000,7) .ALPHA, BETA 
INTEGER  N 


*  PRINT  BANNER  AND  START  SUBROUTINE 
PRINT* 

PRINT*, ******  RUNNING  SIMULATION  LOAD  SUBROUTINE  ****** 
PRINT* 

DO  100  I  =  1,  N 


*  BREAKING  STRENGTH  EQUATION  BASED  ON  WEI  BULL  ALPHA  AND  BETA  PARAMETERS 
DATA( 1,4)  =  EXP( (LOG( -LOG( 1-DATA( 1,3)) )+ALPHA*LOG( BETA) ) /ALPHA) 
100  CONTINUE 
END 
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*  SUBROUTINE  SORT 

Vr  Vr 


SUBROUTINE  SORT  USES  A  BUBBLE  SORT  METHOD  TO  SORT  THE  FILAMENTS 


*  FROM  LOWEST  TO  HIGHEST  BY  BREAKING  POINT  DISPLACEMENT. 

* 


Vr 
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SUBROUTINE  SORT( DATA ,  N ,  E ,  L) 


*  DECLARATION  STATEMENTS 

REAL  DATA(0;  1000,7),T4,T5,E,L 
INTEGER  N, PAIRS, T1 
LOGICAL  DONE 


*  PRINT  BANNER  AND  START  SUBROUTINE 
PRINT* 

PRINT*,’*****  RUNNING  SORT  SUBROUTINE  ****** 

PRINT* 

PAIRS=  N-l 
DONE  =  .FALSE. 

20  IF( . NOT.  DONE)  THEN 
DONE  =  .TRUE. 

DO  30  I  =  1, PAIRS 

IF(DATA( 1,4)  .  GT.  DATA( 1+1,4))  THEN 
T4  =  DATA( 1,4) 

T1  =  DATA( 1,1) 

DATA( 1,4)  =  DATA( 1+1,4) 

DATA( 1,1)  =  DATA( 1+1,1) 

DATA( 1+1,4)  =  T4 
DATA( 1+1,1)  =  T1 
DONE  =  .FALSE. 

ENDIF 

30  CONTINUE 

PAIRS  =  PAIRS  -  1 

GOTO  20 

ENDIF 

PRINT* 

PRINT*,'***  SORT  COMPLETE,  NEW  INDEX  BEING  ADDED  TO  DATA***' 
PRINT* 

DO  40  I  =  1,  N 
DATA( 1,2)  =  I 
40  CONTINUE 
PRINT* 

PRINT*,'***  SOLVING  FOR  DISPLACEMENT  BASED  ON  BREAK  LOADS  ***’ 
PRINT* 

DO  50  I  =  1 ,N 

DATA( 1,5)  =  (DATA( I ,4)*L)/(E) 

50  CONTINUE 
END 
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*  SUBROUTINE  ANLDAT 

Vr  Vr 


*  SUBROUTINE  ANLDAT  WRITES  THE  DISPLACEMENT  AND  LOAD  OF  EACH  BREAK  * 

*  INTO  ARRAY  ADATA 
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SUBROUTINE  ANLDAT( DATA , N ) 


*  DECLARATION  STATEMENTS 
REAL  DATA( 0; 1000,7) 
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INTEGER  N 


*  PRINT  BANNER  AND  START  SUBROUTINE 
PRINT* 

PRINT*, ******  RUNNING  ANALOG  DATA  FILE  SUBROUTINE  ***** ’ 
PRINT* 

DO  100  I  =  1,  N 

DATA( 1,6)  =  (N-I+1)*DATA( 1,4) 

DATA( 1,7)  =  (N-I)*DATA( 1,4) 

100  CONTINUE 
END 


*  SUBROUTINE  DISDAT  * 

*  * 

*  SUBROUTINE  DISDAT  CREATES  LOAD  AND  DISPLACEMENT  BASED  ON  THE 

*  DATA  BUT  THE  DISPLACEMENT  IS  CONTROLLED  BY  THE  MACHINE  CROSSHEA 

*  SPEED  AND  THE  DATA  SAMPLING  RATE.  A  RANDOM  NUMBER  GENERATOR  IS  * 

*  USED  TO  SOLVE  FOR  THE  LOAD  IF  THE  DISPLACEMENT  OCCURS  AT  A  BREAK  * 

*  POINT.  * 

*  * 
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SUBROUTINE  DISDAT( DDATA , DATA , N , E , L , CHS , SR , POINTS ) 
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*  DECLARATION  STATEMENTS 

REAL  DDATA(0;  5000,2) ,DATA(0;  1000,7) ,E, CHS, SR, DELD.DISP, LIMIT, 
1  RANGE, RNG.L 

INTEGER  SEED1, BREAKS, POINTS, N 

*  PRINT  BANNER  AND  START  SUBROUTINE 

PRINT* 

FEINT*,'*****  RUNNING  DISCRETE  DATA  SUBROUTINE  ****** 

PRINT* 

PRINT*.' INPUT  ANY  ODD  INTEGER  (RNG  SEED)' 

READ* , SEED1 
DELD  =  CHS/( SR*60) 

BREAKS  =  0 


*  SET  LIMIT  EQUAL  TO  DISPLACEMET  OF  FIRST  BREAK 
LIMIT  =  DATA(1,5) 

DO  100  I  =  1.5000 
DISP  =  DELD* I 
150  CONTINUE 

IF( DISP  . LT.  LIMIT)  THEN 
DDATA(I.l)  =  DISP 

DDATA( 1,2)  =  ((N  -  BREAKS)*E*DISP)/L 

GOTO  199 

ENDIF 

IF(DISP  .EQ.  LIMIT)  THEN- 

RANGE  =  (((N  -  BREAKS )*E*DI SP) -( ( N -BREAKS - 1 )*E*DISP) ) /L 
DDATA( 1,2)  =  ( RNG ( SEED 1 ) *RANGE ) +( ( ( N - BREAKS - 1 ) *E*D I SP ) / L ) 
DDATA(I.l)  =  DISP 
GOTO  199 
ENDIF 

BREAKS  =  BREAKS  +  1 

IF( BREAKS  .EQ.  N)  GOTO  200 


*  INCREMENT  TO  NEXT  BREAK  DISPLACEMENT 
LIMIT  =  DATA( BREAKS+1 , 5 ) 

GOTO  150 

199  CONTINUE 
100  CONTINUE 

200  CONTINUE 


*  POINTS  =  NUMBER  OF  DISCRETE  DATA  POINTS 

POINTS  =  I 
END 

*  SUBROUTINE  NOISE  * 

*  * 

*  SUBROUTINE  NOISE  ALLOWS  THE  USER  TO  SIMULATE  MACHINE  NOISE  IN  THE  * 

*  DISCRETE  LOAD  DATA.  THE  USER  INPUTS  THE  RANGE  OF  THE  NOISE  AND  * 

*  A  RANDOM  NUMBER  GENERATOR  IS  USED  TO  FIND  A  NOISE  FOR  EACH  DATA  * 

*  LOAD  IN  BETWEEN  THE  GIVEN  RANGE.  * 

Vc  Vc 
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SUBROUTINE  NOISE(DDATA ,  POINTS) 


*  DECLARATION  STATEMENTS 

REAL  DDATA(0; 5000,2), PH, PL.RNG 
INTEGER  I, Q.SEED2, POINTS 


*  PRINT  BANNER  AND  START  SUBROUTINE 
PRINT* 

PRINT*,’*****  RUNNING  NOISE  SUBROUTINE  ***** ' 

PRINT* 

PRINT*,' DO  YOU  WANT  TO  ENTER  NOISE  INTO  DISCRETE  DATA* 

PRINT*,’ IF  YES  TYPE  1,  IF  NO  TYPE  0’ 

READ* , Q 

I F( Q  .NE.  1 ) GOTO  200 

PRINT*,' INPUT  ANY  ODD  INTEGER  (RNG  SEED)' 

PRINT* 

READ*,SEED2 

PRINT*, 'INPUT  HIGH  NOISE  RANGE  AND  LOW  NOISE  RANGE  (GMS) ' 

READ*,  PH,  PL 
DO  100  I  =  0,  POINTS 

DDATA( 1,2)  =  (DDATA( 1,2)  -((PH-PL)/2) )  +  RNG(SEED2)*(PH-PL) 

100  CONTINUE 
200  CONTINUE 
END 

Vr********************************************************************** 


*  SUBROUTINE  OUTPUT  * 

Vc  •>- 

*  SUBROUTINE  OUTPUT  WRITES  THE  ANALOG  DATA  AND  THE  DISCRETE  DATA 

*  ONTO  DRIVE-B  OUTPUT  FILES  * 

Vr  Vr 
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SUBROUTINE  OUTPUT( DATA , DDATA ,N,E,L, ALPHA , BETA , POINTS , ADATA) 
REAL  DATA(0; 1000,7) ,DDATA(0; 5000,2) ,E .ALPHA, BETA, L, 

1ADATA( 0; 2000,2) 

INTEGER  N , POINTS ,Q1 ,Q2 , I , II 

PRINT*,’*****  RUNNING  OUTPUT  SUBROUTINE  ****** 

PRINT* 
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PRINT*, 'DO  YOU  DESIRE  A  LISTING  OF  THE  ANALOG  DATA' 
PRINT*, 'IF  YES  TYPE  1,  IF  NO  TYPE  ZERO' 

READ*,Q1 

PRINT* , ' DO  YOU  DESIRE  A  LISTING  OF  THE  DISCRETE  DATA’ 
PRINT* , ' IF  YES  TYPE  1,  IF  NO  TYPE  ZERO' 

READ*,Q2 

IF(Q1  .EQ.  1)  THEN 

*  WRITE  THE  ANALOG  DATA 
PRINT* 

OPEN(15,  FILE=' ADATA' ) 

WRITE( 15 ,  960) 

960  FORM AT (  /  '  DISP  LOAD’) 

DO  105  I  =  1 ,N 

II  =  ( 2*1 ) -1 

ADATA ( 11,1)  =  DATA(I ,5) 

ADATA( II+l , 1)  =  DATA( 1,5) 

ADATA( 11,2)  =  DATA( 1,6) 

ADATA( 11+1,2)  =  DATA( 1,7) 

105  CONTINUE 

DO  106  I  =  0 , 2*N 

WRITE( 15 , 961)  ADATA(I.l),  ADATA( 1,2) 

106  CONTINUE 

961  FORMAT(  F10. 6,2X,F10. 4) 

CLOSE( 15) 

PRINT* 

PRINT*,'  ANALOG  DATA  IS  IN  FILE  "ADATA  FILE’" 

PRINT* 

END  IF 

IF(Q2  .EQ.  1)  THEN 


*  WRITE  THE  DISCRETE  DATA 

OPEN(15,  FI LE=' DDATA' ) 

WRITE( 15,951) 

951  FORMAT(/  'POINT  DISP  LOAD') 

DO  110  I  =  0, POINTS 

WRITE( 15 ,901)1,  DDATA C 1,1) ,DDATA(I,2) 

110  CONTINUE 

PRINT*,'  DISCRETE  DATA  IS  IN  FILE  "DDATA  FILE’" 
PRINT* 

CLC'3E(  15) 

END  IF 

901  F9RMATC  I5,2X,F10. 6,  4X,F10.4) 

END 
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FUNCTION  RNG(SEED)  RANDOM  NUMBER  GENERATOR 


* 


VcVrVr  “Wr  VrVr  V:  *  *.Wr  VrsWrVr  Vc  V?  ■ 


VpVfVcVrVfVfVcVrVcVc'jV'jVVrVf'jV^V^VVc'jVVrVr'sVVr*’ 


Vr*  Vr  Vr  VnWr  VoVtV^V-jWc*  VrVc  Vr^VVr  Vf -sWr  iV^WcVriV 


FUNCTION  RNG(K) 

INTEGER  K,M, CONST 1 
REAL  RNG ,  CONST2 

PARAMETER(CONSTl=2147483647 ,  CONST2=. 4656613E-9) 
SAVE 

DATA  M  /0/ 

IF(M  .EQ.  0)  M  =  K 
M  =  M*65539 

IF(M  .LT.  0)  M  =  (M+l)  +  CONST 1 
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RNG  =  M  *  C0NST2 
END 
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APPENDIX  C 


COMPLIANCE  REMOVAL 

The  objective  of  the  analysis  of  fiber  bundle  testing  data  is  to 
determine  the  fiber  strength  distribution  parameters.  To  this  end.  the 
actual  fiber  displacement  data  must  be  obtained  from  the  total  test 
displacement  data.  A  significant  portion  of  the  test  displacement  data 
is  due  to  machine  compliance.  This  appendix  describes  a  method  to 
determine  this  displacement  and  subtract  it  from  the  test  displace¬ 
ment  data  consisting  of  the  load  train  displacement  and  the  sample 
displacement. 

A.  ISOLATION  OF  COMPLIANCE 

Fiber  bundle  testing  is  accomplished  with  the  use  of  the  Instron 
Universal  Testing  Instrument  Model  4206  (INSTRON)  machine.  The 
fiber  test  records  displacement  and  the  load  measured  by  the 
INSTRON  load  cell.  The  test  is  controlled  by  software  that  operates 
the  machine  on  a  controlled  cross- head  speed  and  takes  data  at  a 
user-selected  sampling  rate.  The  machine  is  schematically  illustrated 
in  Figure  43. 

A  typical  bundle  test  load  vs  displacement  graph  is  shown  in  Fig¬ 
ure  44.  The  curve  may  be  partitioned  into  three  regions.  Region  I  is 
the  concave  non-linear  region.  This  is  concavity  is  due  to  the  compli¬ 
ance  on  the  load  train  and  the  variation  in  gauge  length  (slack)  amount 
of  each  filament. 
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Figure  44 

Typical  Bundle  Test  Load  vs.  Displacement  Graph 
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This  is  due  to  each  filament  of  the  bundle  not  having  exactly  the 
same  gauge  length;  the  curve  becomes  progressively  steeper  as  each 
filament  starts  to  carry  the  load.  Region  II  is  the  linear  region.  Here  all 
of  the  filaments  are  carrying  the  load  equally.  The  bundle  displace¬ 
ment  is  a  linear  function  of  displacement  due  to  Hooke’s  Linear 
Stress-Strain  law  shown  in  Equation  23. 

P  =  nE5 

5  =  P/nE  (23) 


where  P  =  load 

n  =  number  of  filaments 
E  =  filament  modulus 
5  =  displacement 

Region  III  is  the  failure  region  where  the  load  fluctuates  as  each 
fiber  breaks. 

The  load  data  that  the  INSTRON  records  is  derived  from  an  elec¬ 
trical  signal  from  the  load  cell.  This  load  is  the  load  taken  by  the 
fibers,  but  also  includes  the  load  taken  up  by  the  testing  machine  load 
train.  This  load  causes  deformations  in  the  load  train  consisting  of  the 
gears  of  the  machine,  the  grips,  grip  base,  and  all  other  structural 
components.  These  components  of  the  load  train  combined  to  form 
machine  compliance.  The  INSTRON  shown  in  Figure  43  may  now  be 
modeled  as  shown  in  Figure  45.  The  INSTRON  is  equivalent  to  two 
springs  in  series,  the  machine  compliance  and  the  fibers. 
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Figure  45 

INSTRON  Testing  Machine  Equivalent  Model 


With  the  spring  in  series  model,  the  total  displacement  of  the 
system  due  to  a  load  is  equal  to  Equation  24.  The  load  due  to  the  fibers 
is  Equation  25  and  the  load  due  to  compliance  is  equation  26. 


8t  =  Sf+8c  =  P/(Kf+Kc) 


(24) 


8f  =  P/Kf 

(25) 

8C  =  P/Kc 

(26) 

where  8t  =  total  bundle  displacement 
8f  =  fiber  displacement 
8C  =  compliance  displacement 
Kf  =  fiber  stiffness  coefficient 
Kc  =  compliance  stiffness  coefficient 
Equations  25  and  26  demonstrate  that  the  fiber  and  compliance 
displacements  are  functions  of  load.  If  the  stiffness  coefficients  and 


96 


the  load  are  known,  then  the  displacement  due  to  the  fibers  may  be 
found  by  Equation  27.  It  is  important  to  note  that  the  stiffness  coeffi¬ 
cients  for  the  fibers  and  the  compliance  are  likely  to  be  nonlinear.  The 
fiber  nonlinearity  is  primarily  due  to  the  slack  in  the  early  displace¬ 
ment.  The  compliance  will  be  nonlinear  due  to  the  physical  processes 
that  contribute  to  form  the  compliance.  Many  of  these  processes  are 
directional  in  nature. 


5(P)f  =  8(P)t  -  8(P)c  (27) 

From  Equation  27,  the  typical  bundle  test  graph  In  Figure  1  may 
now  be  represented  as  Figure  46  with  the  compliance  displacement 
being  subtracted  from  the  fiber  displacement. 


Figure  46 

Typical  Bundle  Test  Load  vs  Displacement  Graph  Compliance 
and  Fiber  Displacement  Separated 
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To  determine  the  fiber  displacement,  a  means  to  subtract  the 
compliance  displacement  is  required.  An  equation  must  be  deter¬ 
mined  to  fit  8(P)C  to  implement  Equation  27. 

If  Kc  were  linear,  it  would  be  a  simple  matter  to  solve  for  the  8(P)C 
function.  However,  as  mentioned  above,  the  function  will  not  be  linear. 
Because  of  this,  8(P)C  must  be  measured  by  experimentation.  The  data 
must  then  be  curve-fit  to  form  a  continuous  function;  this  function  may 
then  be  used  to  implement  Equation  27. 

B.  COMPLIANCE  TESTING 

From  the  model  Illustrated  in  Figure  45,  a  means  to  eliminate  the 
displacement  of  the  fibers  would  be  to  set  the  gauge  length  of  the 
bundle  to  zero.  The  fibers  may  not  be  eliminated  because  the  gripping 
mechanism  to  hold  the  bundle  in  place  is  a  part  of  the  compliance,  but 
the  gauge  length  can  be  brought  very  nearly  to  zero.  To  do  this,  a  sam¬ 
ple  is  prepared  like  the  one  shown  in  Figure  47.  If  there  are  no  fibers 
that  displace,  all  of  the  displacement  recorded  must  be  due  to 
machine  compliance. 

To  measure  8(P)C,  a  bundle  test  is  completed  with  the  0  gauge 
length  sample  in  accordance  with  the  procedures  outlined  in 
Appendix  D.  A  typical  INSTRON  compliance  test  load  vs.  displacement 
graph  is  shown  as  Figure  48. 
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Figure  47 

0  Gauge  Length  Sample  Used  to  Test  Compliance 


Figure  48 

INSTRON  Compliance  Test  Load  vs.  Displacement  Graph 
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The  first  region,  0  to  Del,  is  very  nonlinear.  This  is  due  to  the 
directional  compliance  in  the  INSTRON  as  the  displacement  is  first 
applied.  Physical  processes  include  the  backlash  in  the  gears,  nonlin¬ 
ear  friction  and  damping  as  the  gears  and  clutches  engage,  and  noise 
below  the  sensitivity  of  the  load  cell. 

Region  two,  Del  to  Dc2,  is  the  linear  region.  The  compliance  in 
this  region  is  due  to  linear  spring  tension  of  the  machine  as  modeled 
in  Figure  44.  The  K  is  very  high,  which  is  expected  for  a  testing 
machine.  It  should  be  high  so  that  the  load  it  absorbs  is  very  low  when 
compared  to  the  load  being  absorbed  by  the  sample. 

Region  three  is  nonlinear.  It  is  believed  that  in  this  region  the 
grips  are  beginning  to  fail  and  that  the  fibers  are  actually  slipping 
between  the  copper  tabs.  The  data  in  this  region  does  not  reflect  the 
machine  compliance  and  is  not  utilized  in  the  analysis.  It  is  assumed, 
however,  that  if  the  sample  had  not  slipped  the  compliance  curve 
would  have  continued,  extending  the  linear  region. 

C.  CURVE-FITTING  PROCEDURE 

The  fitting  of  a  curve  to  approximate  P(6)c  was  completed  by 
separating  the  function  into  two  regions,  selecting  a  power  series 
model.  Equation  28,  for  region  I,  a  linear  model;  Equation  29.  for  sec¬ 
tion  II;  and  setting  boundary  conditions  at  Del,  Equations  30  and  31. 


for  0  <  Dc  <  Del 

Dc  =  Al(P)n  +  A2 

(28) 

for  Del  <  Dc 

Dc  =  A3P  +  A4 

(29) 
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Boundary  Conditions 


(1)  Del  (Equation  28)  =  Del  (Equation  29)  (30) 

(2)  Del'  (Equation  28)  =  Del'  (Equation  29)  (31) 

If  Dc2  and  Del  are  known,  then  from  Equation  29  for  region  two: 

A3  =  Dc2  -  Dcl/Pi  -  p(Dcl)  (32) 

A4  =  -A3P1  +  Del  (33) 


From  Equations  28,  29,  and  31,  an  equation  for  Ai  can  be  derived: 

8  Dc/8  P  =  nAi(P)n-l 
8  Dc/8  P  *  A3 
A3  =  nAifP)0'1 

A2  =  A3/n  (PiJn-l  (34) 

From  Equations  29,  30,  and  34,  an  equation  for  A2  can  be  derived: 

Del  =  ((A3/n  (Del  -  DcO)n-i)  *  (Del  -  Dc0)«/A3  +  A2 

A2  =  Dei  -  A3Pi/n  (35) 

This  gives  us  five  unknowns  and  four  equations.  Equations  32  and 
33  can  be  solved  directly  from  information  gained  from  the  0  gauge 
length  data.  Equations  34  and  35  Eire  coupled  by  n,  so  A3  and  A4  may 
be  solved  for  by  optimizing  the  variable  n.  To  optimize  n,  an  error 


function  is  created  and  its  derivative  is  taken  and  then  summed  for  all 
the  data  points  between  zero  and  Del  for  different  values  of  n.  When 
the  error  function  is  minimized  to  approximately  zero,  n  is  known. 
The  n  can  then  be  substituted  into  equations  13  and  14  to  yield  the 
constants  A3  and  A4. 

Substituting  Equations  32  and  33  into  Equation  28  gives  Equation 
36. 

Dc(P)  =  (A3Pi/n)  *  ((P/Pi)«-l)  +  Del  (36) 

3  Dc(Pi)/3  N  =  (A3Pi/n)  *  ln(P/Pi)  *  (P/Pi)n 

+  A3Pin-2(l  -  (P/P!)n)  (37) 

The  error  function  becomes 

H(n)  =  X  (Dc(Pj)  -  DO  3  Dc(Pi)/a  N  (38) 

for  i  =  1,2,3  ...  number  of  data  points  from  0  to  Del.  Equation  38  is 
then  iterated  until  H(n)  goes  to  zero  (.0000). 

The  process  of  solving  for  the  coefficients  is  greatly  enhanced  by 
the  use  of  spreadsheets  and  microcomputers.  Programs  such  as  Lotus 
1-2-3  and  Microsoft  Excel  allow  the  user  to  import  the  data  directly 
into  a  spreadsheet  and  manipulate  the  data.  The  results  for  the 
INSTRON  compliance  used  in  this  project  were  found  using  Lotus 
1-2-3. 
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A  Fortran  program  is  in  development  to  solve  this  procedure 
is  not  functional  at  this  time. 


APPENDIX  D 


BUNDLE  TESTING 

The  failure  load  for  graphite  bundles  was  determined  with  the  use 
of  the  Instron  Universal  Testing  Instrument  Model  4206  (INSTRON) 
and  associated  IBM  INSTRON  control  software,  version  4.01.  The  fiber 
strength  statistics  were  found  by  removing  the  compliance  displace¬ 
ment  with  the  software  developed  in  this  appendix.  The  compliance 
displacement  function  was  approximated  by  the  curve  fit  detailed  in 
Appendix  C.  The  modified  data  was  then  used  with  software  developed 
by  Lt.  Joseph  Schmidt  to  determine  the  fiber  parameters  a  and  p  for 
the  given  gauge  length. 

All  testing  was  completed  in  the  Advanced  Composites  Laboratory, 
Naval  Postgraduate  School,  Monterey.  Laboratory  Technician  Mr.  Jim 
Nageotte  and  Professor  Edward  Wu  assisted  with  the  preparation  of 
samples  and  the  testing. 

A.  BUNDLE  TESTING  PROCEDURE 

The  testing  was  completed  in  accordance  with  the  following 
procedure. 

•  Turn  on  INSTRON  and  allow  60  minutes  to  warm  up. 

•  Calibrate  the  INSTRON  in  accordance  with  posted  procedures  for 
mechanical  and  electrical  calibration  (only  required  for  first  test 
of  the  day). 

•  Turn  on  attached  IBM  AT  and  enter  INSTRON  test  software.  From 
C  prompt,  type:  MT\DATA\MT  <enter> 

•  Load  sample  into  machine. 
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Ensure  sample  is  vertically  aligned  correctly;  adjust  if  necessary. 
Set  zero  load  level. 

Depress  Load  Bal  button  and  then  enter  on  the  INSTRON. 

Set  cross-head  speed  to  .2  mm/min. 

Set  gauge  length. 

Run  cross-head  up  until  a  load  of  .2  ±.02  kg  is  achieved. 
Depress  Gauge  Length  button  on  the  INSTRON. 

Run  cross-head  down  until  a  load  of  zero  is  achieved. 

Run  cross-head  up  until  a  load  of  .2  ±.02  kg  is  achieved. 
Extension  should  read  0.0.  If  it  does  not,  repeat  all  steps. 

Run  cross-head  down  until  a  load  of  zero  is  achieved. 

Depress  jog  down  switch  twice. 

Depress  Gauge  Length  button. 

"(Crosshead  Speed  may  be  adjusted  if  necessary) 

Select  appropriate  test  method. 

Enable  IEEE  interface. 

Depress  IEEE  button  on  the  INSTRON. 

Run  the  test. 

Plot  the  results. 

Plot  on  HP  plotter  using  Plot  option  in  INSTRON  software. 

Put  data  into  ASCII  file. 

Dump  data  to  hard  disk  or  floppy  disc  using  Utilities  option  in 
INSTRON  software. 

Exit  INSTRON  software. 

Remove  Compliance  displacement. 

Rename  test  data  file  to  TEST.DAT  and  put  onto  floppy  disk. 

Place  data  disk  in  drive  B: 
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•  Place  compliance  removal  program  disk  in  drive  A: 

•  Type  T5  <enter> 

The  disk  in  drive  B:  now  contains  the  data  with  compliance 
removed  and  the  data  with  compliance  not  removed.  The  data  is  now 
ready  to  be  run  in  the  software  developed  by  Lt.  Joseph  Schmidt  that 
will  retrieve  the  fiber  parameters  a  and  (5. 

B.  PROGRAM  LISTING 

The  compliance  removal  program  “program  T5"  takes  the  five 
coefficients  for  the  curve  fit  compliance  function  and  then  for  the  cor¬ 
responding  load  subtracts  the  displacement  from  the  data  due  to 
compliance.  The  method  for  determining  the  compliance  curve  fit 
coefficients  is  detailed  in  Appendix  C. 

The  program  is  written  in  Microsoft  Fortran  4.01.  It  will  run  on 
any  IBM/compatible  personal  computer.  To  operate  the  program, 
place  the  program  disk  into  the  A:  drive.  The  data  disk  must  be  placed 
into  the  B:  drive.  When  this  is  completed,  type:  T5  <enter>. 

PROGRAM  T5 

SUBROUTINE  LISTING: 

I  NIT  Initializes  arrays  to  zero 

LOAD  Opens  the  data  file  and  reads  in  the  machine  data 

CONVERT  Converts  machine  data  into  SI  units 

REMCON  •  Prompts  the  user  to  input  the  compliance  curve 

coefficients 

•  Removes  the  displacement  from  the  data  due  to 
compliance 


OUTPUT 


Writes  to  the  B:  drive  the  modified  data 


U- 


VARIABLE  LISTING: 

DDISP(6000) 

Test  displacement  data 

DLOAD(6000 

Test  load  data 

INDEX 

Test  data  point  counter 

XCONV 

Displacement  conversion  factor 

YCONV 

Load  conversion  factor 

A1 

Compliance  curve  coefficient 

A2 

Compliance  curve  coefficient 

N 

Compliance  curve  coefficient 

A3 

Compliance  curve  coefficient 

A4 

Compliance  curve  coefficient 

DC1 

Transition  displacement  between  two  curve  fit 
equations 

ZERO 

Noise  level  of  load  cell 

n 


Bj 
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PROGRAM  T5 

***Vr**Vr*Vr********7V*^V*Vc********'i^*^V*Vc**iV**************Vr******'5V***5V**** 

*  PROGRAM  T5 

*  * 

*  THIS  PROGRAM  WILL  TAKE  INSTRON  4200  SERIES  BUNDLE  RAW  DATA  OUTPUT  * 

*  AND  ALLOW  THE  USER  TO  CONVERT  THE  DATA  INTO  DESIRED  UNITS  AND  TO  * 

*  TO  REMOVE  MACHINE  COMPLIANCE  DISPLACEMENT.  THE  INPUT  DATA  FILE  * 

*  MUST  BE  PLACED  IN  THE  B-  DRIVE.  THE  INPUT  DATA  FILE  MUST  BE  * 

*  NAMED  TEST. DAT.  A  MAXIMUM  OF  6000  DATA  POINTS  IS  ALLOWED.  TO  * 

*  EXECUTE  THE  PROGRAM,  PLACE  THE  PROGRAM  DISK  IN  THE  A  DRIVE,  * 

*  TYPE  A-T5  AND  HIT  ENTER.  tHE  USER  SHOULD  THEN  FOLLOW  THE  SCREEN  * 

*  PROMPTS.  * 

*  THE  REQUIRED  INPUTS  ARE-  * 

*  -  NUMBER  OF  DATA  POINTS  * 

*  -  CONVERSION  FACTOR  FOR  LOAD  AND  DISPLACEMENT  * 

*  -  COMPLIANCE  CURVE  PARAMATERS  (IF  DESIRED)  * 

*  * 

*  THE  OUTPUT  DATA  FILE  WILL  BE  PUT  ON  THE  DISK  IN  THE  B  DRIVE  AS 

*  FILE  EXPER. OUT.  A  SUMMARY  PRINTOUT  OF  THE  INPUT  IS  PLACED  ON 

*  THE  DISK  IN  THE  B  DRIVE  AS  DATCON.  PRT  .  * 

*  * 

REAL  XCONV,YCONV,Al,A2, A3, A4,N,DDISP( 6000), DLOADC 6000), 

+ZERO 

INTEGER  I , J,DP, INDEX( 6000) ,Q 
PRINT*,'***  RUNNING  PROGRAM  T5  ***' 

PRINT* 

*  CALL  SUBROUTINE  INIT  TO  INITIALIZE  ARRAYS 

CALL  INIT( DDISP, DLOAD, INDEX, XCONV, YCONV) 

*  CALL  SUBROUTINE  LOAD  TO  LOAD  IN  DATA  INTO  PROGRAM 

CALL  LOAD( DDISP , DLOAD , INDEX . DP ) 

*  ALLOW  USER  TO  CONVERT  DATA  INTO  DESIRED  UNITS 

PRINT*,'  DU  YOU  NEED  TO  CONVERT  DATA  INTO  [SI]  UNITS' 

PRINT*,'  THIS  SHOULD  ONLY  BE  REQUIRED  FOR  MACHINE  TEST  DATA' 
PRINT*,'  IF  YES  TYPE  1,  ELSE  TYPE  0  TO  CONTINUE' 

READ* , Q 

IF(Q  .EQ.  1)  THEN 

*  IF  USER  DESIRES  TO  CONVERT  UNITS  CALL  SUBROUTINE  CONVERT  TO  CONVERT 

*  DATA 

CALL  CONVERTC XCONV ,YCONV, DP, DDISP, DLOAD ) 

END  IF 

*  CALL  SUBROUTINE  REMCOM  TO  REMOVE  COMPLIANCE  FROM  DATA 

CALL  REMCOMC DDISP, DLOAD , A1 , A2 , A3 , A4 ,N , DP , ZERO , DC  1 ) 

*  CALL  SUBROUTINE  OUTPUT  TO  OUTPUT  DATA  INTO  DATA  FILES  ON  DISK  B 

CALL  OUTPUT ( INDEX , DD I SP , DLOAD , A 1 , A2 , A3 , A4 , N , DP , ZERO , XCONV , YCONV , 
+DC1) 

PRINT*,'***  END  OF  PROGRAM  ***' 

END 

*********************************************************************** 

*  SUBROUTINE  INIT  * 

*  * 

*  THIS  PROGRAM  WILL  INITIALIZE  THE  DATA  ARRAYS,  INDEX  AND  CONVERSION  * 

*  FACTORS  PRIOR  TO  IMPUTING  NEW  DATA. 

*  is 

SUBROUTINE  INIT( DDISP, DLOAD, INDEX, XCONV, YCONV) 
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REAL  DDISP( 6000)  ,DLOAD(6000) , XCON'V , YCONV 
INTEGER  INDEX( 6000) 

*  PRINT  HEADER 

PRINT* ,  '***  RUNNING  SUBROUTINE  INIT  ***' 

PRINT* 

*  RUN  LOOP  TO  SET  DATA  ARRAY  INDEX,  DDISP  AND  DLOAD  TO  0 

DO  10  I  =  1,6000 
INDEX(I)  =  0 
DDISP(I)  =  0. 

DLOAD(I)  =  0. 

10  CONTINUE 

*  SET  CONVERSION  FACTORS  TO  1.  0 

XCONV  =  1. 

YCONV  =  1. 

END 

*********************************************************************** 


*  SUBROUTINE  LOAD 

* 


* 


*  SUBROUTINE  LOAD  WILL  OPEN  DATA  FILE  B-TEST.  DAT  AND  READ  IN  THE  * 

*  DATA  WITH  AN  UNFORMATTED  READ  STATEMENT.  THE  NUMBER  OF  DATA  * 

*  POINTS  (DP)  TO  BE  READ  MUST  BE  ENTERED  CORRECTLY  OR  PROGRAM  WILL  * 

*  ERROR  OUT  IN  THIS  SUBROUTINE. 


Vr 


* 

* 


SUBROUTINE  LOAD( DDISP , DLOAD , INDEX , DP ) 

REAL  DDISP(6000) ,DL0AD(6000) 

INTEGER  I ,DP , INDEX(6000) 

PRINT  HEADER 

PRINT*,’***  RUNNING  SUBROUTINE  LOAD  ***’ 

PRINT* 

OPEN  DEVICE  15,  FILE  B-TEST. DAT 
OPEN( 15 , FI LE=’ B-TEST.  DAT' ) 

INPUT  NUMBER  OF  DATA  POINTS  (DP) 

PRINT*, ' INPUT  NUMBER  OF  DATA  ELEMENTS  TO  BE  READ’ 
PRINT*,’  MAXIMUM  INPUT  VALUE  IS  6000  (2(5(3 ’ 

READ* , DP 

RUN  LOOP  TO  READ  IN  DATA  POINTS  INDEX,  DDISP  AND  LOAD 
DO  10  I  =  1 ,DP 

RE AD ( 1 5 , * ) I NDEX( I ) , DD I SP ( I ) , DLOAD ( I ) 

10  CONTINUE 
CLOSE  DEVICE  15 
CLOSE  (15) 

END 
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*  SUBROUTINE  CONVERT  * 

*  * 

*  SUBROUTINE  CONVERT  WILL  CONVERT  THE  DATA  FROM  FILE  TEST.  DAT  INTO  * 

*  WHICHEVER  UNITS  THE  USER  DESIRES.  THE  USER  MUST  INPUT  CONVERSION  * 

*  FACTORS  WHICH  WILL  BE  DIVIDED  FROM  THE  ORIGINAL  DATA.  THE  X  CONV  * 

*  CORRESPONDS  TO  DISPLACEMENT  AND  THE  Y  CONV  CORRESPONDS  TO  LOAD  OR  * 

*  FORCE.  IF  NO  CONVERSION  IS  DESIRED,  ENTER  1  FOR  THE  CONVERSION 

*  FACTORS. 


SUBROUTINE  CONVERT( XCONV , YCONV ,DP , DDISP , DLOAD) 
REAL  DDISP(6000) ,DL0AD( 6000) , XCONV, YCONV 
INTEGER  I, DP 
*  PRINT  HEADER 
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PRINT*,'***  RUNNING  SUBROUTINE  CONV  ***' 

PRINT* 

PRINT*,  'PROGRAM  T5  WILL  CONVERT  DATA  INTO  [SI]  UNITS' 

PRINT* 

*  INPUT  X  OR  DISPLACEMENT  CONVERSION  FACTOR 

PRINT* , ' INPUT  X  CONVERSION  (DISPLACEMENT)' 

READ*,  XCONV 

*  INPUT  Y  OR  LOAD/FORCE  CONVERSION  FACTOR 

PRINT*, 'INPUT  Y  CONVERSION  (LOAD/FORCE) ' 

READ*,  YCONV 

*  CONVERT  DATA 

DO  10  I  =  1 ,DP 
DDISP(I)  =  DDISP( I) /XCONV 
DLOAD(I)  =  DLOAD( I ) /YCONV 
10  CONTINUE 
END 

*Vf*Vf*****Vf*************Vr*****,iWfVpiWriV****Vf  ******************************* 

*  SUBROUTINE  REMCON  * 

*  * 


*  SUBROUTINE  REMCON  WILL  IF  THE  USER  DESIRES  REMOVE  THE  MACHINE  * 

*  COMPLIANCE  DISPLACEMENT  FROM  THE  DATA.  THE  USER  MUST  INPUT  THE 

*  COMPLIANCE  CURVE  PARAMETERS  AND  THE  DATA  ZERO  LOAD  LEVEL. 

*  THE  USER  SHOULD  TAKE  CARE  TO  ENSURE  THE  COMPLIANCE  PARAMETERS  ARE  * 

*  THE  SAME  UNITS  AS  THE  CONVERSION  FACTORS. 


SUBROUTI NE  REMCOMC DD I SP , DLOAD , A 1 , A2 , A3 , A4 , N , DP , ZERO , DC 1 ) 

REAL  DDISP( 6000) ,DLOAD(6000) ,A1 , A2 , A3 , A4 ,N , ZERO ,DC1 
INTEGER  I ,DP,Q 

*  PRINT  HEADER 

PRINT*,'***  RUNNING  SUBROUTINE  REMCOM  ***' 

PRINT* 

*  ALLOW  USER  TO  REMOVE  COMPLIANCE  FROM  THE  DATA 

PRINT*, 'DO  YOU  DESIRE  TO  REMOVE  THE  MACHINE  COMPLIANCE  FROM  THE  DA 
+TA' 

PRINT*, 'IF  YES  TYPE  1,  ELSE  TYPE  0  TO  CONTINUE' 

READ*,Q 

IF(Q  .EQ.  1)  THEN 

*  IF  USER  DESIRES  TO  REMOVE  COMPLIANCE,  INPUT  COMPLIANCE  PARAMETERS 

PRINT*, ' INPUT  COMPLIANCE  CURVE  COEFFICIENTS' 

PRINT*, 'INPUT  Al’ 

READ*, A 1 

PRINT*, ' INPUT  A2 ' 

READ* , A 2 
PRINT* , ' INPUT  N' 

READ*,N 

PRINT*, ’ INPUT  A3' 

READ*, A3 

PRINT*, ’ INPUT  A4 ' 

READ* , A4 

PRINT* , ' INPUT  DC1 ' 

READ*, DC 1 

PRINT*, ' INPUT  ZERO  NOISE  LEVEL  (LOAD  IN  KG)' 

READ*, ZERO 

*  REMOVE  COMPLIANCE  DISPLACEMENT  FROM  THE  DATA 

PRINT* 

PRINT*,'***  REMOVING  MACHINE  COMPLIANCE  FROM  DATA  ***' 


l  in 


PRINT’’' 

DO  100  I  =  1  ,DP 

*  DETERMINE  IF  DISPLACEMENT  DATA  IS  IN  COMPLIANCE  POWER  CURVE  OR 

*  LINEAR  REGION  AND  SUBRTRACT  CORRECT  AMOUNT  OF  DISPLACEMENT 

IF(DLOAD(I)  .  LT.  ZERO)  THEN 
GOTO  99 


END  IF 

IF(DDISP(I)  .LT.  DC1) 
DDISPd)  =  DDISP(I)  - 


THEN 

((Al’v(DLOADCI)**N))  +  A2) 

ELSE 

DDISPd)  =  DDISP(I)  -  ((A3’vDLOAD(I))  +  A4) 

END  IF 

99  CONTINUE 
100  CONTINUE 
END  IF 
CONTINUE 
END 

*  SUBROUTINE  OUTPUT  * 


*  SUBROUTINE  OUTPUT  WILL  WRITE  THE  DATA  TO  FILE  EXPER. OUT.  IT  WILL  * 

*  ALSO  CREATE  A  DATA  SUMMARY  FILE  AS  DATCON. PRT  .  THE  OUTPUT  FILES  * 

*  ARE  PLACED  ONTO  THE  DISK  IN  THE  B  DRIVE.  * 

V«V*VfV**VVfVfifiV),VVfVrVfV?VrVfVfVfVfiVVf>VVrVfVfVfVfVfVfVfV:VfyfVfVfVrVfVf»Vf^VfVfVoV:VVnVVfTV*Vf?'c:VVoV*V>V,j,fVc«VVfVfVc^**,fVf»Vyf^ 


SUBROUTINE  OUTPUT( INDEX , DD I SP , DLOAD , A 1 , A2 , A3 , A4 , N , DP , ZERO , XCONV , 
+YC0NV.DC1) 

REAL  DDISPC 6000) ,DL0AD( 6000) ,A1 ,A2 , A3 ,A4 ,N , ZERO , XCONV , YCONV 
INTEGER  I,DP,INDEX(6000),Q 

*  PRINT  HEADER 

PRINT*,’***  RUNNING  SUBROUTINE  OUTPUT  ***' 

PRINT* 

10  PRINT*, 'DATA  WILL  BF.  OUTPUT  TO  FILE  B -EXPER.  OUT’ 

*  PROMPT  USER  TO  PLACE  TARGET  DATA  DISK  IN  DRIVE  B- 

PRINT*, ’ENSURE  A  DISK  WITH  SUFFICIENT  SPACE  IS  IN  DRIVE  B-’ 
PRINT*, 'TYPE  1  TO  CONTINUE' 

READ* , Q 

IF(Q  . NE. 1)  GOTO  10 

*  OPEN  DEVICE  12  AND  WRITE  DATA  TO  FILE  B -EXPER.  OUT 

OPEN( 12 ,  FILE=' B -EXPER. OUT’ ) 

DO  100  I  =  1 ,DP 

WRITE( 12, 900) INDEXC I), DDISF(I ),DLOAD(I) 

100  CONTINUE 

*  CLOSE  DEVICE  12 

CLOSE( 12) 

*  ALLOW  USER  TO  CREATE  A  DATA  SUMMARY  FILE 

PRINT*, 'DO  YOU  DESIRE  A  SUMMARY  PRINTOUT  FOR  THE  DATA  CONVERSION* 
PRINT*, 'IF  YES  TYPE  1,  ELSE  TYPE  0  TO  EXIT’ 

READ* ,Q 

*  IF  USER  DESIRES  DATA  SUMMARY  FILE,  OPEN  DEVICE  12  AND  WRITE 

*  PARAMETERS  TO  FILE  DATCON. PRT 

IF( Q  . EQ.  1)  THEN 

OPEN(  12 ,  FILE**'  B -DATCON.  PRT'  ) 

WRITE( 12,910) 

WRITE(12,911)DP, XCONV, YCONV 
WRITE ( 12 ,912)ZER0,DC1 
WRITE( 12 ,913)A1 ,A2 ,N,A3,A4 


CLOSE( 12) 

END  IF 
CONTINUE 

*  WRITE  FORMAT  STATEMENTS 

900  FORMAT(  16, IX, F9. 4, 1X.F9. 4) 

910  FORMAT(  'COMPLIANCE  REMOVAL  DATA  CONVERSION  COEFFICIENTS',//) 

911  FORMAT (  'NUMBER  OF  DATA  POINTS' ,8X, 1 7/ 

1  'X  CONVERSION  FACTOR  ' ,5X,F10. 7/ 

2  'Y  CONVERSION  FACTOR  ' ,5X,F10. 7//) 

912  FORMAT (  'ZERO  NOISE  LEVEL  (KG) ' , 10X.F7. 4/ 

1  ' DC1  -LINEAR  SECTION-  '.8X.F9.4//) 

913  FORMAT(  'COMPLIANCE  CURVE  COEFFICIENTS’// 

1  'POWER  FACTOR  PORTION'// 

2  ' A1  =  ' ,F6. 4/ 

3  ' A2  =  ' ,F6. 4/ 

4  '  N  =  ,F6. 4/ 

5  'LINEAR  PORTION'// 

6  'A3  =  ,F6. 4/ 

7  'A4  =  ' ,F6. 4/) 

END 


APPENDIX  E 
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TEST  DATA 

A.  COMPLIANCE  TEST  DATA 

The  INSTRON  compliance  test  data  forms  are  contained  in  this 
appendix.  They  are  tabulated  in  Test  Number  sequence.  The  raw  data 
is  stored  on  the  IBM  AT  personal  computer  in  the  Composites  Labora¬ 
tory,  Naval  Postgraduate  School,  in  director  MT/DATA/  as  file  (Test 
Number). MAD. 


J 
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NFS  Compcsites  Lab 

U.S. Naval  Postgraduate  School 

Monterey,  CA 


Fiber  Bundle  Compliance  Curve 
2ero  guage  length 
.15  aa/am 

Test  type:  Compliance  Curve  I  ns  t  r  on  Corporation 

Spmps  II  Automated  Materials  testing  System  v4.A|c 
Operator  na*e:  Jim  Test  Date:  August  5,19E9 


Sample  Identif ication:  0GL1C5 
Interface  Type:  4200  Series 
Machine  Parameters  of  test: 

08 

Sample  Type: 

bndl 

Saaple  Rate  (pts/sec): 

10.00 

Haeidity  1  « 

1:  60 

Crcsshead  Speed  (mi/min 

):  .150 

Temperature  (deg.  F): 

90 

Di  /rensions: 

Spec.  1 

Lin.  Density  (der) 

5036 . 0 

Gauge  length  < e* ) 

i .c;:o 

Cut  of  1  specimens,  0  excluded. 

Lead 

Displcrent 

I  Strain 

Displceenl 

lead 

hctfalus 

at 

at 

at 

at 

at 

Spec i een  label1 

Pariepe 

Brut 

Brea! 

Maximum 

Ka»ir-.» 

Nueter 

(ge/denl 

(Kg) 

Ire) 

m 

(Ml 

(Kg) 

1  gcc  d 

97.95 

29.57 

.3011 

30.11 

.£942 

£9.59 

Mean: 

97.95 

£9.57 

.3011 

30.11 

.2992 

£5.59 

Standard 

Deviation: 


Mean  - 
2.00  »  Sdv: 


Mean  ♦ 

2.00  *  5dv: 


Minimum: 

C3.95 

£9.57 

.3011 

3C.11 

•  £9.£ 

2"  .57 

Man  mu*: 

9?.  95 

£9.57 

.3011 

3. \ll 

,£9i£ 
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vuuiiiuN  a .  »:>  j  r 


Bon.  5 1 

1  (-’  1  (1 

OGL 1 ODOO 

Tr*-.t  date 

August  V>,  1900 

Vpr r.  j  on 

A  .Ole 

Version  elate 

06/07/00 

Rcpor t 

f  1  1  t?tt 

3A 

Operate r 

J  j  m 

X  convnr j on 

.  <03937*  >o& 

X  A/D  ofO.Pt 

.  (  H'l(n'l 

V  convt?r*»icn 

P- .  20A62300 

Y  A/D  of! set 

.  *  >ooo 

Samp  1  e 

r  a  t  c 

1 0 . 00 

E>:  t  rnsc 'meter 

NO 

A/D 

range? 

o 

Autc*r.  tart 

or  r 

Cal  ib 

t  ypr 

AU10MA1 IC 

t  r  y 

CVL  INDfx  I  (  A! 

CY-  1  i  h 

1  P.ul 

1 10. 2 31 1 

Ca  )  ih  e  :  t*M  i«\ 

.  <.»'  * 

T  r-mpc*r 

r*  t  ure 

70 

Humid) t  y 

60 

i^t 

type? 

or;  n  i  t;r< 

tt  specimens 

i 

Din- 

type/ 

bnd  1 

Entry  cl  i  mens 

YDS 

E<  r  e>  a  \ 

check 

1  00 .  00000 

Thresh  delay 

.0 

toil  cl 

l  i  m  i  t 

1 10.231 10 

l  ::  t  enr>  1  l  n.  i  t 

.  i : -voo 

Sump  1  e 

d  i  mei  «p  i  .on**  : 

A: 

.  039 A 

D  :  2036 .  0000 

C: 

. 039A  D: 

. 0390  E :  NO 

Bp’ or  ] 

mr-'n  0 

i 

Step  roll,  status  :  PKt.AI  Di  li’ 

Mo  ::  )  (Tii.ii! 

)  e*(1 

60 . 191 

Ma::  I  cad  point  tl  :  1  i  T<0 

Mo;:  & 

tt-.-.s 

.013 

Max  extent  pint  :  1309 

Numbrr  r>  f  c  le»«ent  £•  :  1309 

Bpr.i.  )  »iif ')  i  cJ  i  iim:  j  on*-,  : 

A  :  .  039370 1  l * :  3030 . 0'  *00'.»o0  CJ  :  .  037370  ID:  .  03V 37v )  T  :  N1 ) 


Att::j  1  i  ;  ry  Input  Array  1 

**<*»*<  <«**♦•»*  ♦  .*■*■«**¥  •*#»»**  «  *«***<-«  ««*«-««•  ##««««*  «  y  »  «  •  «  N  *  *  M  *  »  «  * 

Air.'  i  liar  y  Input  Army  3 

*  *  *  ■«  4  *  *  «***.**..*  *+■¥****  *-•****♦.  #**-*-««*-  *•»*»#*#  *  M  «  «  «  *  *  *«*■»»*■»  .  i;  I  ( 
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B.  BUNDLE  TEST  DATA 

The  INSTRON  bundle  test  data  forms  are  contained  in  this 
appendix.  They  are  tabulated  in  Test  Number  sequence.  The  raw  data 
is  stored  on  the  IBM  AT  personal  computer  in  the  Composites  Labora¬ 
tory.  Naval  Postgraduate  School,  in  directory  MT\DATA\  as  file  (Test 
Number). MAD. 


CPMFOS 1 TES  LftBOfiAlORY 
NAVAL  POSTGRADUATE  SCHOOL 
HONTEREY  CA  S3940 

GRAPHITE  BUNDLE  TEST  HETH30  (70  CHI 
Ust  aettiod  for  instron  4306.  Uied  for  testing  AS-4 
Graphite  bundle!  Ijppi  1000  fiber!) 

Test  type:  OIL  BUNDLE 

Operator  naoe:  HARK  JONES 

Sanple  Identification:  080901 
Interface  type:  4300  Seriei 
Nachine  Parameter!  of  test: 

Saaple  Rate  I pts/sec ) :  6.67 

Crosehead  Speed  lai/am  ):  7.000 


‘ BATCH  ID'  060901  70  CN  WITH  OIL 

Dieeneion!: 

Spec.  I 

Lir.  Deriity  Iter)  5036.0 

6au;.  .enyL.  (at)  700.00 

Out  of  1  speciaens,  0  excluded. 

Displcaent  J  Strain 
at  at 

Speciaen  Naxiaua  Naxiaua 
N'jaber  leal  (t) 


1 

5.240 

.71.86 

.135? 

Hejn: 

5.240 

.7486 

.1359 

Sta"dj'd 

Deviaticn: 


Load 

at 

Naxiaua 

If.N) 


(nation  Corporation 

Sene!  11  Autoaated  Hater  i a  1  s  Testing  Systea  v4.0lc 
Test  Date;  Septeaber  8,1939 

Saaple  type:  AS-N 

Huaidity  III:  SO 

Teaperature  Ideo.  F 1 :  73 
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COMPOSITES  LABORATORY 
naval  postgraduate  school 

MONTEREY  CA  939NO 


GRAPHITE  BUNDLE  TEST  15  CHI 


Test  type:  Yarn/fiber 

Operator  na»e:  Mark  Jones 


Instron  Corporation 

Series  II  Automated  Materials  Testing  Systee  vN.Olc 
Test  Date:  September  9,1988 


Sa*ple  Identification:  090901  Saiple  Type:  AS -4 

Interface  Type:  4200  Senes 
Machine  Parameters  of  test: 

5i»plf  Rate  (pts/sed:  3.33  Huaidity  (  J):  50 

Crosshead  Speed  (u/ein  h  .500  Teeperature  (deq.  FI:  73 


'BATCH  ID'  090701 

'OIL  OR  NO  Oil1  NO  OIL 

Dimensions: 

Spec.(l 

Lin.  Density  (teil  2034.0 

6auge  length  (•«)  50.000 

Out  of  1  specimens,  0  er.tluded. 


Displceent 

X  Strain 

Lead 

at 

at 

at 

Spec i»pn 

Mauiv* 

Han  pup 

Nanny* 

N-jrber 

(ma> 

(1) 

(KN) 

1 

.6500 

!  .300 

.2349 

Mean; 

,6!"0 

1.300 

.2349 

Standard 

Deviation: 


fl^an  - 
2.00  i  Sdv: 

Mean  * 

2.00  i  Sd<: 


Minieu*; 

.6500 

1.300 

.2349 

Ma*ieu»: 

.6500 

1.300 

.2349 
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VERSION  4.01c 


Sample?  id 

070 VO 1 

Test  date 

September  7,  lVfUl 

Vf?r  ni  on 

6  .  I.l  1  c 

Version  date 

06/09/00 

F<epcTt  f  ilc  » 

12 

Operator 

Mori-  .To nrs 

X  convprr.ion 

. 03937000 

X  A/n  off st?t 

.  00IW| 

V  cr«nveraif«n 

E  .  20'*  62300 

Y  A/I)  of  f  set 

.  !.)«.)»  >'  ) 

Samp  It**  ralp 

3.33 

Ex  tonsomo ter 

NO 

A/D  r  a  nci  e 

o 

Autostar  t 

OFF 

Cal  lb  typo 

AUTOMATIC 

Geometry 

CYLINDRICAL 

C  a  1  )  h  1  o  n  d 

1 10.231 1 

Ca  1  j  h  ex  tpnr» 

.  0000 

Tempera  lure? 

73 

Hum i d i t  y 

SO 

Test  type 

DFNIER 

♦l  specimens 

1 

F<ar  type 

AS -4 

Entry  d l mens 

YES 

Eire;d  clir-nl 

1 0 . 0000*.) 

Thresh  delay 

.  P.t'UV  1 

l.  o  a  d  1  i  m  \  t 

1 12. A 0630 

Ex  tens  limit 

.  OVll'i  .1 

Sample  dimenc 

ions  : 

A :  .  0394 

D:  10326.0000  C: 

1 . 9605  D: 

.0396  F  i  1 10 

Spec). mm  <1 

1 

Test  end  status  :  10 

Maximum  load 

02.703 

Max  load  point:  1)  :  .'•/»*> 

Max  o.xtcns 

.047 

M a  x  extens  pnt  H  :  4 A 3 

Number  of  elements  :  4C>3 

Spot-  i  men  d  j  mi'Mir  1  cn  .r.  : 

As  . 0373701  1«:  103:?'* .  00o0«  )«.»0  C:  1 . 960304 «.»  ))  :  .Ol-JVPVnj  (  :  Nil 


Auxiliary  Input  Arr  */  tt  1 

«4**»**  -**««.«#«  *#*.*♦*.¥  *  *  ■*  ■<  #  **  *  **<#.*#  * « # « » * *  *•**♦**«  #.###*■«*  ******* 

Auxiliary  Input  ttrrr-y  <1  P 

#  ^*  *  *  *  *  #***■-***  *******  *  «-««**  W  *«■*■**•**  *******  *******  *******  **»*«*« 
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COMrOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939^0 


GRAPHITE  BUNDLE  TEST  (.5  Cm 


Test  type:  Yarn/Fiber 

Operator  name:  Hark  Jones 


Irstron  Corporation 

Series  IX  Automated  Materials  Testing  System  v 4 . 0 1 c 
Test  Date:  September  9, 1903 


Sa.-nle  Identification:  090902 
Interface  Type:  4200  Series 
Machine  Parameters  of  test: 

Sample  Rate  (pts/sec):  .33 

Crosshead  Speed  (ma/tin  ):  .050 


Sample  Type:  AS-4 

Humidity  (  I  ):  50 

Temperature  (deg.  F):  73 


•BATCH  ID'  090902 

m  OR  MOIL'  NO  OIL 


Dimensions: 

Spec.  1 

Lin.  Density  (ter)  2C2E.O 

Gauge  length  (am)  5.0000 

Out  of  I  specimens,  0  excluded. 


Oisplcmert  X  Strain  Load 
at  at  at 

Specimen  Maximum  Maximus  Maximum 
Number  (ms)  IX)  .  fkN) 


I  .2500  5.000  .2034 


Mean:  .2500  5.000  .2034 

Standard 

Deviation:  -  -  - 


Mean  - 
2.00  *  Sdv: 

Mean  ♦ 

2.00  *  Sdv: 


Minimum: 

.2500 

5.000 

.  203s* 

Maximum: 

.2500 

5.000 

.2034 
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COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93940 


GRAPHITE  BUNDLE  TEST  (.5  CM) 


Test  type:  Yarn/P » ber 

Operator  nare:  Mark  Jones 


Instron  Corporation 

Series  IT  Automated  Materials  Testing  System  v4.01c 
Test  Date:  September  9 ^ J 983 


Sample  Identification:  090903 
Interface  Type:  4200  Series 
Machine  Parameters  of  test: 

Sample  Rate  Ipts/sec):  .33 

Crosshead  Speed  (••/■in  ):  .050 


Sarple  Type:  AS-4 

Huiidity  (  J  ):  50 

Teiperature  (dpg.  F):  73 


'PATCH  ID'  090903 

'OIL  OR  NO  OIL1  NC  Oil 

Dimensions: 

Spec.  1 

Lin.  Density  (tpx)  2036.0 

Gauge  length  («*)  5.0000 

Out  of  1  specimens,  0  ercluded. 

Sarple  co>rents:  Mann*  Load  22.36  Tg  Relation  Load  to  21.33  Yq 


Oisplc^ent 

X  Strain 

,  Load 

at 

at 

at 

Specimen  Ma*i»ut 

Ma*ieu* 

Harisui 

Number  (••) 

(I) 

ON) 

1  .2500 

5. 000 

,?w 

Mean:  ,2500  5.000  .2194 

Standard 

Deviation:  . 


Mean  - 
2.00  »  Sdv: 

Mean  ♦ 

2.00  *  Sdv: 


Nmi»u*: 

.2500 

5.000 

Man ; *u* : 

.?5r0 

5.000 

.2W 

I 
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COMPOSITES  LftSOPflTORY 
KflVflL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93940 


GRAPHITE  BUNDLE  TEST  (.5  CHI 


Test  t ,-pe:  Yarn/Fiber 

Operator  name:  Hart  Jones 


Instron  Corporation 

Series  II  Automated  Materials  Testing  Systee  v4 . 0 1c 
Test  Date:  September  9,1988 


Saeple  Identification:  090909 
Interface  Type:  8800  Series 
Machine  Parameters  of  test: 

Saeple  Rate  (pts/sec):  .33 

Crosshead  Speed  lpn/*in  ):  .050 


Saeple  Type:  AS-8 

Humidity  III:  50 

Temperature  ideg.  FIs  73 


'BATCH  ID'  090908 

'OIL  OR  NO  OIL'  NO  OIL 

Dimensions: 

Spec.  1 

Lin.  Density  (tei)  8036.0 

Gauge  length  Imml  5.0000 

Out  of  I  specimens,  0  eicluded. 

Sample  comments:  Specnin  is  slipping  and  test  stepped  at  .96  u 


Di  sp  1  client 

S  Strain 

trad 

at 

at 

at 

Specimen 

Hu  taut 

Hatitum 

Mantu* 

Nu*ber 

(•ml 

m 

t  ►  N ) 

I 

.8700 

5.800 

.8357 

Mean: 

.8700 

5.800 

.8357 

Standard 
Devi  at  ion: 


Mean  - 
8.00  »  Sdv : 

Mean  ♦ 

8.00  e  Sdv: 


Him  tun: 

.8700 

5.800 

.2357 

Han*ut: 

.8700 

5.800 

.2357 
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COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  939^0 


GRAPHITE  BUNDLE  TEST  1.5  CM) 


Test  type:  Yarn/Fiber 

Operator  na*t:  Nark  Jones 


Instrpn  Corporat ion 

Series  II  Autonateo’  Materials  Testing  System  .01 
Test  Date:  September  9, 19E9 


Sample  Identification:  090905 
Interface  Type:  4300  Series 
Machine  Parameters  of  test: 

5a»ple  Rate  tpts/sec);  .33 
Crosshead  Speed  (••/•in  ):  .050 


Saiple  Type:  AS-4 

Humidity  t  X  I:  50 

Temperature  (deg.  FI:  73 


'PATCH  ID'  090905 

’OIL  OR  NO  OIL’  NO  OIL 


Dnensicns: 

Spec.  1 

Lin.  Density  (te*)  3036.0 

6auge  length  (n)  5.0000 


Out  &f  1  specimens,  0  excluded. 
Sample  ccseents:  Saeple  slipped 


D  i  sp lc»ent 

X  Strain 

.  Itid 

at 

at 

at 

Sf?ci«e<! 

Hariri* 

Maritue 

Max  i  run 

N-jrbfr 

tea) 

(X) 

UNI 

1 

.8500 

5.800 

.8998 

Mean: 

Standard 

Cevuticn: 

.8900 

5.800 

.8998 

Mean  - 

3.00  *  Sd v ; 

Mean  ♦ 

3.00  *  Sdv: 

Mini«ue: 

.3900 

5.800 

.3993 

Max i rut: 

.3900 

5.800 

.3992 
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COMPOSITES  LABORATORY 
NAVAL  P0ST6R4DUATE  SCHOOL 
HONTEREY  CA  939«.0 


GRAPHITE  BUNDLE  TEST  (.5  CHI 


lest  type:  Yarn/Eiber 

Operator  ease:  Hark  Jones 


I  ns  t  r  on  Corporation 

Series  II  Automated  Hjleruls  Testing  System  vk.Olc 
Test  Date:  September  9,1988 


Slfple  Identification:  090904 
Interface  Type:  4309  Series 
Hacnme  Pirneters  of  test: 

Satple  Rate  Ipts/secl:  .33 

Crossbead  Speed  (me/min  I:  .050 


5a«ple  Type:  AS-4 

Hjeidity  I  t  I:  50 

Tearerature  ( deq .  FIs  93 


'BATCH  10'  09.191,4 

■OIL  OR  NO  OIL'  HO  OIL 

Dimensions: 

Spec.  1 

Lin.  Density  (teil  3034.0 

Gauge  length  5.0000 

Out  of  1  sptcieens,  0  evduded. 

Satole  toeients:  SAMPLE  SLIPPED 

t  Strain  .  load 

at  at 

Haiimum  AaiieuA 
III  (PHI 


Oisplctent 

at 

Srecioen  Hane-jo 
Nuiber  (os) 


1  .3400  A. BOO  .3444 


Hean:  .3400  4.800  .8444 

Standard 

Deviation:  -  -  - 


bear  - 

?. 00  *  Sdv: 

Hein  ♦ 

3.00  I  Sdv: 


Hiniaua: 

.2400 

4.  BOO 

.3444 

Pjiiaua: 

.2*00 

4.B00 

.3444 
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COMPOSITES  LABORATORY 
NAVAL  F0ST6^ACUAT£  SCHOOL 
HONTEREY  Cfl  93940 


GRAPHITE  BUNDLE  TEST  (5  CM 


Test  type:  Yam/Fiber 

Opf'ator  na*e:  Hark  Jones 


Instron  Corporation 

Series  II  Automated  Materials  Testing  System  v4.pjc 
Test  Date:  Septeeber  9,j9E3 


Saeple  Identif nation:  090907  Saeple  Type:  A3-4 

Interface  Type:  4800  Series 
Machine  Para*eters  of  tost: 

Sample  Rate  ipts/secl:  3.33  H;*idity  i  s  j:  cn 

Crosshead  Speed  ta«/ain  )•  .500  Treperature  f deg .  F!:  73 


'PATCH  IP'  090707 

•OIL  Cfi  NO  OIL'  NO  OIL 

[l»Pr3|C,nt; 

Spec.  I 

Lip.  Density  lie*)  3'“ 36.0 

6 i j 9 e  Tenet's  <m»  50.000 


Out  cf 
Saiple 


1  spec A 
CCf*?nt$:  S^“:lE 


excluded. 

SLIFFED  AT  END  O'  TEST  TEST 

Cl  Spicier.' 

at 

Specimen  Manege 
NtmTer  leal 


STORED  F*  C-fEFATCP 

t  X  St*  am  .  lead 

at  at 

Muciaua  Ha*  i  *•„•* 

(1)  (M.n 


1  .6500  1.300  .20  36 


Mean:  .fcCOv  1.300  .Vet 

Stag'd 

Deviation:  - - - 


H«=an  - 

2.0*'  •  Sd v : 

Hpj"  ♦ 
p.on  a  3d v : 


Hi  n  i  : 

.6500 

1.300 

.2026 

Ma-ieut: 

.  blCO 

1.3'Vi 

.2026 

i 


125 


CCHPOSilES  LSPC'SAIORY 
HfiVPL  fOSISSACUATE  SCHOOL 
HONTEREY  CA  939AO 


GRAPHITE  BUNDLE  TEST  (5  CHI 


I  ns  tron  Corporation 

Senes  II  Autoeated  Materials  looting  Systr.  v 
lost  Catei  Septesfcer  9,1985 

Saaple  Typo:  AS-A 


Hjeidity  (  t  I:  50 

Teeperature  (deg.  FI:  73 

'PATCH  IS'  090908 

'OIL  08  NO  OIL'  NO  OIL 

DlPSnsions: 

Spec.  t 

Li”.  Density  <  t * 1 1  2035 . 0 

taug;  length  (eel  50.000 

Ojt  of  I  epecieer.,  A  erdoded. 


lest  type:  Yarn/Fiber 

Operator  pane:  Hart  Jones 

Saeple  Identification:  09 09na 
Interface  Type:  9300  Series 
Hachine  Faraeeters  of  test: 

Satple  Rate  (ptsfsed:  3.33 
Crosshead  Speed  (ee'nn  I: 


Di srl c»ert 

J  Strain 

lead 

at 

at 

at 

Waxi*ji 

Navi  tut 

Nusber 

(«M 

m 

( KN  > 

) 

.fc4(  0 

1.280 

.3:99 

Wea^: 

.6400 

1 .250 

.3339 

Standard 

Deviation: 


Hean  - 

3.00  *  Sdv ; 

3.00  i  Sd>: 


Piniui: 

.6400 

1.380 

.3393 

Wax 

.6400 

1.350 

.3333 
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COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93940 


GRAPHITE  BUNDLE  TEST  15  CM! 


Test  type:  Yarn/Fiber 

Operator  r,a?e:  Mark  Jones 


Jnstron  Corporation 

Spries  II  Automated  Materials  Testing  System  v4.01c 
Test  Cate:  September  9,1989 


Sample  Identification:  090909 
Interface  Type:  4800  Series 
Machire  Paraeeters  of  test: 

Sample  Rate  (pts/sec):  3.33 
Crosshead  Speed  (ii/im  I:  .500 


Sample  Type:  AS  4 


Hueidity  IX):  50 

Teeperature  (deg.  F):  73 


’BATCH  It’  090909 

'OIL  OR  NO  OIL*  OIL 


Dimensions: 

Spec.  1 

Lin.  Density  flex)  8036.0 

Gauge  length  (eei  50.000 


Out  cf  I  specimens,  0  preluded. 
Sample  cos'ents:  OIL  SAMPLE 


Displcuert  X  Strain  }.rad 
at  at  at 

Sp?cu?n  Mausum  fla*if -ji  Mane.* 

N-'ter  1m!  ( r. )  (MI) 


1  .6500  1.300  .8374 

Mean:  .6500  1.300  .8374 

Standard 

Deviation:  . - . . . 


Mean  - 
8.00  *  5d v : 

Mean  ♦ 

2.00  *  Sdv: 


Ninitut: 

.6500 

1.300 

.337'. 

ManU'jt: 

.65r'0 

1.300 

.3374 
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VERSION  4.01c 


Samp  )  e  id 

090909 

T  c  s  t  date 

5r*ptpmlmr  9,1720 

Version 

A  .  0  1  c 

Version  date 

06/09/88 

Report  filed! 

12 

□per  a to r 

Marl  Jonr& 

X  conversion 

. 03937008 

X  A/D  offset 

.  nt  ii  ii  > 

Y  conversion 

2.20A62300 

Y  A/D  offset 

.  I."  >0*  ) 

Samp  If’  i'A  1 1? 

3.33 

E::  tensomotor 

NO 

A/D  range 

0 

Autostart 

OFF 

Ca lib  type 

AUTOMATIC 

Geomo  try 

CYLINDRICAL 

C  a  ]  i  b  load 

110.231 1 

Ca  1  ib  ex  tens 

,  000o 

1  PinpnrA  t  i.tr  e 

73 

Hum  i  rj  i  ty 

50 

Test  type 

DCNI Eft 

#  spec i mens 

1 

E«ar  type 

AS-  A 

Entry  d i mens 

YES 

Breal  chrcl 

1  0  .  OOt  >00 

Thresh  delay 

.  PEA  FT  1 

Load  limit 

112.  AT)  A  CO 

E ::  tens  limit 

.09  f  i-i  3 

Sample  fJjnir*nsic*i>5 
A:  .039'*  E* 


Spec  i  n i c • » •.  If 
M • '  i  mi  im  1  ( ■  | 

1*1. V 


1 83c?4 . 0000  C: 


J 

• ,:i .  pvt, 
.  <.>SV 


i  .  9605  I): 


"If?  ?»t  end  r.tatu 
M<a::  lead  point 
M  a  :<  p : :  t  f  .m  i  o  pul 


.0390  FI:  r  in 


-  :  1  «.» 

II  :  .  ;v,i 

U  :  4*  ’J 


Nuff.hr r  of  t-  l  ♦  men  l  s  :  603 


Sp« 'cinipn  dimensions  : 

A:  .0393701  h:  1  H3P4 . 0('»0000O  C  :  1.96US040  b:  .0393  701  Fit!'’ 


Auxiliary  Input  Army  II  1 

*«*¥«««  «•«*«»«  ¥•**¥•***  ¥¥«***#  «+««*.««  *¥»«««¥  ¥¥*♦*.«« 

Auxiliary  Input  Array  II 

¥**¥**¥  *¥¥»*¥*  *¥»*¥¥¥  *«««*«#  MM«M 


««««■«¥  ¥*¥**#* 
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COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
HQNTEREt  CA  <53940 


GRAPHITE  BUNDLE  TEST  (2.5  CHI 


Test  type:  YarnfFiber 

Operator  na«e:  Hark  Jcnes 


I  ns  t  r  on  Corporation 

Series  I)  Automated  Materials  Testing  Systee  v4.0lc 
Test  Oate:  September  9,1988 


Sanple  Identification;  090910 
Interface  Type:  4200  Series 
Machine  Parameters  of  test: 

Saeple  Rate  (pts/sec > :  1.47 

Crosshead  Speed  (mm/min  I:  .250 


Saeple  Type:  flS-4 

H'joidity  It);  50 
Teeperature  (deg.  E):  73 


'BATCH  10'  070510 

’OIL  OR  NO  OIL'  NO  OIL 


Diffusions: 

Spec.  1 

Lin.  Density  (te»>  2034.0 

Gauge  length  lee)  25.000 

Out  of  1  specimens,  0  excluded. 


Displcfsnt 

X  Strain 

Lead 

at 

at 

,  ct 

Specif?n 

Pantus 

Naxivuc 

Haiirjt 

N'jrber 

(m) 

(X) 

UN) 

1 

.4600 

1 .840 

.2179 

^5  an : 

.4600 

1.350 

.5174 

Standard 

Deviation: 


Mean  - 
5.00  •  Sdv: 

Mean  ♦ 

5.00  f  Sdv: 


Riniauc: 

.4600 

1.840 

.5174 

Hautus: 

.4600 

t  .840 

.2174 
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COMPOSITES  LPBOPAIOPY 
NAVAL  POSIGBflCL'AIF  SCHOOL 
MONTEREY  CA  93940 


GRAPHITE  BUNDLE  TEST  (2.5  CHI 


Test  type:  Yarn/Fiber 

Operator  nj«e:  Hark  Jones 


lnstron  Corporal  ion 

Series  II  Automated  Materials  Testing  Systea  vt.Olc 
Test  Date:  Septeaber  10,1968 


Saeple  Identification:  100901 
Interface  Type:  4200  Series 
Machine  Parameters  of  test: 

Sample  Rate  (pts/secl:  1.47 
Crosshead  Speed  (aa/ain  ):  .250 


Sample  Type:  AS-4 

Huaidity  I  I  I:  50 

Temperature  (deg.  FI:  73 


'BATCH  ID'  100901 

'OIL  OR  NO  OIL'  NO  OIL 

Oieensions: 

Spec.  I 

Lm.  Density  Ite.i  2026.0 

Gauge  length  (tel  25.070 

Out  of  1  «pec i sens ,  o  excluded. 


Displcaent 

X  Strain 

Lc-ad 

at 

at 

at 

5p?cir?n 

Kaxicui 

Haxiaua 

Maxima 

Nj«tpr 

(■*) 

(X) 

(MO 

1 

.4300 

1.720 

.2130 

tfean: 

.4350 

1.720 

.2130 

Standard 

Deviation: 


Kean  - 
2.00  <  Sdr: 

Mean  * 


Mmrja: 

.4300 

1.720 

.2130 

llaxieut : 

.4300 

1.720 

.2130 
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x>  cn 


VERS ION  4 . O 1 c 


Sample*  id 

100901 

7  e s t  date 

September  10, 1 7f 

Vers i on 

4 . 0 1  c 

Version  da'.c? 

06/09/82 

Report  file# 

14 

Opes- a  t  or 

Mas-!  Jones 

X  conversion 

. 03737000 

X  A/D  of  f  4jf' t 

.  0<  »UM 

Y  conversion 

B. 80463300 

Y  A/D  of  1 5ft 

,  O'  s'.jO 

Oaf  up  l  r*  rate 

1  . 67 

E ::  t  mrsDniP  i.  r»t 

NO 

A  / 1)  ra  ng  e 

o 

Au  t  c*  s  t  a  r  t 

cirr 

Ca lib  type 

AUTOMATIC 

Geomi?  t  r  y 

CYLINDRICAL 

Ca 1 ib  1 oad 

110.2311 

Calib  e>: tens 

.  0000 

Tempera  ture 

73 

Hum  l  d  i  1;  y 

50 

Tr-r^t  type? 

DENIER 

specimens 

1 

Bar  type 

AS -4 

En  t  r  y  d  i  mes  is 

YES 

Brfial:  chock 

1 0 -  00000 

Thresh  delay 

1  .  12404 

Lear)  limit 

.1  12 .40450 

E::  tens  limit 

.<14921 

Sample  dimenn 

ions  : 

A :  . 0394 

B  :  1 2324 . 0000 

C: 

.9043  D: 

.0394  E :  NO 

Spec i men  ft 

1 

Test  end  slati 

is  :  1  <> 

Mn  \  rrtt  sm  1  r<  nd 

4  7 .  n:  v  > 

M  ■' ::  l  i."Hl  po  i  n 

II  :  in'! 

Matt  c* t. «  ..s 

.  Oi  '7 

Nat:  f::tc*u«.  pn 

,  tl  : 

Number  of  elements  :  305 

S  p  tv  c  i  ns  e  n  dimension «-  : 

A  :  .  03737V 1  f.< :  1  83r?4 .  vWnVon  f.  :  .  704  f LT-Eu  D  :  .  039'*r?0  1  f  t  Mf  i 


Austi  1  lary  Input  As  ray  t»  1 

*  -#  <  *  *  V  *  »***¥**•  -4  »¥  *  **  *  *  «  «  «  K  *  »  *******  *  *  *  *  *  *  *  *******  *******  ******* 

Au::i  1  i  o?ry  Input  Array  K  P 

*******  M  *M«IM  *¥**.**#  *******  ¥  ¥  *.  *  *  *  *  *******  ******* 
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COflfOSTTES  laboratory 
NAVAL  POSTGRADUATE  SCHOOL 
HON1ERET  CA  73760 


GRAPHITE  BUNDLE  TEST  (2.5  CHI 


Test  type:  Yarn/Fiber 

Operator  name:  Hark  Jones 


1  ns t r on  Corporation 

Series  II  Automated  Haterials  Testing  Systee  vA.Olc 
Test  Date:  September  10,1 98B 


Saeple  Identification:  1 00702 
Interface  Type:  R 200  Series 
Hachine  Paraeeters  of  test; 

Sasple  Rate  ( pts/sec ) :  1.67 

Crosshead  Speed  (»e/«in  ):  . 250 


Saeple  Type:  AS-A 

Hundity  (II:  SO 

Temperature  (deg.  Eli  73 


'PATCH  ID'  1 0070! 

■OIL  OR  NO  OIL'  OIL 

Dimensions: 

Spec.  I 

Lin.  Density  (te»l  5036.0 

Gauge  length  (is)  55.000 

Out  of  1  specimens,  0  ercluded. 

Saeple  comments:  SAMPLE  APPEARED  TO  SLIP 


Di5plc»ent 

X  Strain 

Load 

at 

at 

at 

$ppcii»en 

flam'ja 

Haxiauf 

Hanifj* 

NuutJer 

(Ml 

(X) 

r  (IN) 

i 

.5000 

5.000 

.1755 

Hean: 

.5000 

5.000 

.1755 

Standard 

Deviation: 


Neae  - 
5.00  •  Sdv: 


MpVI  ♦ 

2.0-1  •  Sdv:  . . . 


Kmia'j-: 

■  50C0 

e.cco 

.1755 

Harieut: 

.50f,0 

2.000 

.1755 
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COUPOSITES  LABORATORY 
naval  rosigraouaie  school 
honterey  CA  93940 


GRAPHITE  BUNDLE  TEST  (2.5  HI) 


Test  type:  Yarn/Fiber 

Operator  na«e:  flark  Jones 


Instron  Corporation 

5eries  II  Animated  tlaterials  Test i ng  Systee  vt.Olc 
Test  Date:  Septeeber  10, I960 


Saeple  Identification:  100903 
Interface  Type:  4200  Series 
Hactiine  Paraoeters  of  test: 

Saeple  Rate  (pts/sec):  1.47 
Crosshead  SpeeJ  (et.'mn  ):  .250 


Saople  Type:  AS-4 

Hoaidity  (lit  50 
Teaperature  (beg.  Tl:  73 


'BATCH  ID'  103903 

'OIL  OR  NO  OIL'  OIL 

Cie.ensions: 


Spec.  I 


Lin.  Density  (te>l  2034.0 

Gauge  length  (to)  25.000 


Out  of  I  spec  leers,  0  eicloOei. 
Sa'ple  coseents:  SMflE  SLIPPED 


Dispic*?rt 

X  Strain 

Lead 

at 

at 

•  at 

S;?ci«en 

flarifui 

Narieua 

N‘j«t-er 

(rn) 

(X) 

IFNI 

1 

.4300 

1.720 

.2051 

Itean: 

.4300 

1.720 

.2051 

Standard 

Deviation: 


M 


Nean  - 
2.00  •  Sdv: 


Hojin  * 


8.(0  *  Sdv 


HinjR'j*: 

.4300 

1.720 

.2051 

.4300 

1 .720 

.2051 
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COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93940 


GRAPHITE  BUNDLE  TEST  150  CM) 


Test  type:  Yarn/fiber 

Operator  name;  Mark  Jones 


Instron  Corporation 

Series  II  Automated  Material?  Testing  System  vA.Olc 
Test  Date:  September  10,1 998 


Sample  Identification:  100904 
interface  Type:  4200  Series 
Machine  Parameters  of  test: 

Sample  Rate  (pts/sec):  6.67 
Crosshead  Speed  tmm/min  ):  5.000 


Sample  Type:  A5-4 


Humidity  (  1  >:  50 

Temperature  (deg.  F):  73 


•BATCH  ID'  100904 

'Oil  OR  NO  Oil'  NO  Oil 

Dimensions: 

Spec.  1 

Lin.  Density  ftp*)  2036.0 

Gauge  length  (**)  500.00 

Out  of  l  specters,  0  excluded. 


Di*pk»ent 

1  Strain 

lead 

at 

at 

at 

Specimen 

fla.H'j* 

Maximum 

Maximum 

Number 

(«il 

(1) 

UN) 

i 

3.370 

.6740 

.1070 

Mean: 

3.370 

.6740 

.1070 

Standard 

Deviation: 


Mean  - 
2.00  «  Sdv: 

Mean  ♦ 

2.00  •  Sdv: 


Minimum: 

3.370 

.6740 

,1070 

Maximum: 

3.370 

.67A0 

.1090 
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VERSION  4.01c 


Sample  id 

1 00904 

Test  date 

September  1 u , j 90S 

V trsi o n 

4.01c 

Version  date 

06/09/1)8 

Report  filrft 

1  ' 

Operator  : 

Marl:  .Tones 

X  conversiDn 

*  03737000 

X 

A/D  offset 

V  conversion 

2 . 20462300 

Y 

A/D  offset 

.  O'j'.iO 

Sample  rate 

6.67 

E 

:  tense-meter 

NO 

A/D  range 

0 

Autos tar t 

of  r 

Ca  1  i b  type 

AUTOMATIC 

Geometry 

CYLINDRICAL 

Cal ib  load 

1 10.231 1 

Calib  ertens 

.  000- . 

1 empnr a ture 

73 

Humidity 

50 

Tc?r»t  type 

DEN I E R 

It  specimens 

i 

Par  type 

AS-4 

Entry  dimens 

YES 

Dr  eat.  chect' 

1 0 . 00000 

Thresh  delay 

.22401 

Load  limit 

112. 40450 

E 

:  tens  limit 

.984;".; 

Sample  dimensions  : 

A :  .  0394 

B :  18324.0000 

C: 

19 

.  6850  D : 

.039*4  h:  ML) 

f!p«r.  ininn  II 

1 

Test  end  stati 

•  *  :  1" 

Mn  ::  i  muni  load 

24 . 472 

11 

■> ::  1  ( • .  u  1  p  o  I  1 1 

U  :  r*  7 1 

Ha:: 

.  295 

Ha:;  f.- ::  tens  pri 

t  «  :  6") 

Number  of  e i emmi  t s  :  60 1 

S u  c-  c  i  m r.» n  H  i  nir n s  j  o  n o  : 

A  :  .  0393701  H  :  1  0324  .  OOOOOOO  r. :  19. 6G50400  V  :  .  03937»  >1  f;  :  NM 


Auxiliary  Input  Army  It  1 

***«>*■••  *»*«**«  *  X  *  *  *  *  >«  *  +  *  H  *  +  * 

Auxiliary  Input  Array  It  P 

*  «  V  *  ¥  •*  «  If- *  «  0- #  X  *  *  *  *  *  *  «  »  #  «  «  »  « 


***.*.  X  *  *  «  4  •  *  *  «  X  *  *  4  *  *  >  ►  *  k  *  *  .  *  4 

«  «  V  4  X  4  4  *  *  *  4  4  »  4  X  »  «  «  *  4  »  M»»‘n  *»»»"*» 
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COMPOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
HONTEREY  CA  939AO 


GRAPHITE  BUNDLE  TEST  150  CHI 


Test  type:  Yarn/Fiber 

Operator  naae:  Hark  Jones 


Instron  Corporation 

Series  |(  Autoaated  Haterials  Testing  Systea  vA.Olc 
Test  Date:  September  10,1988 


Saaple  Identification:  100905 
Interface  Type:  ASCO  Series 
Nacbme  Paraeeters  of  test: 

Sa»ple  Rate  (pts/sec):  4.67 
Crosshead  Speed  (••/•in  I:  5,000 


Saaple  Type:  AS-A 


H'jaidity  IT):  50 
Teaperature  Ideg.  FI:  73 


'PATCH  ID  inn?05 

'OIL  OR  NO  OIL'  OIL 

D i Tensions: 

Spec.  I 

Lin.  Density  Iter)  8036.(1 

Gauge  length  (eel  507.00 

Out  of  I  scecieens,  0  excluded. 

Saeple  cements:  OIL  SAftFLE 


Displt»ent 

*  Strain 

Icatf 

at 

at 

at 

Specimen 

hauiui 

Ramus 

Number 

(it) 

(X) 

Ik  N) 

1 

3.910 

.7850 

.139? 

fean: 

3.910 

.7850 

.1399 

Stand-X'll 
Ceviat ion: 


Hean  - 
5.00  •  Sdv: 

Hean  * 

3.00  *  Sdv: 


Mimiu*: 

3.910 

.7830 

.1399 

Na>i«u«: 

3.910 

.7830 

.139? 
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Tost  d?tp 
Version  date? 
Operate*  r 


iirptrmtirr  IVOR 

06/09/00 

Marl  Jones 


VERSION  4.01c 


Sample  id 

100905 

Ver s i on 

4.0lc 

Report  file*! 

1  1 

X  conversion 

. 03937008 

Y  conversion 

2 . 20462300 

Sample  rate 

6.67 

A/D  range 

0 

Cal i b  type 

AUTOMATIC 

Calih  load 

110. S3 1 1 

T  emper a ture 

73 

Test  typo 

DENIER 

par  type 

A3 -4 

break  check 

1 0 . 0000O 

Load  limit 

112. 40450 

Samplp  (1 1  mrnr  i  on'.  : 

A:  .0394 

L» ;  3  S3 P 4 . 0000 

Spec,  j  men  It 

i 

Maximum  It -ad 

31  .47) 

Max  ox  tons 

.314 

Number  of  elements  :  644 

Spec i men  d i mens i ons  : 

A  :  .  03'7CJ7h  lb:  1  0334  .  Om(>  m  u‘»0  C  : 


X  A/D  offset  : 

.  0<  •<><> 

Y  A/D  offset  : 

.  OuOO 

E ::  t  pnsc-neter 

NO 

Autos tar t 

OFF 

Gcoinetr  y 

CYLINDRICAL 

Ca  1  l  b  ox  tons 

.  oooc* 

Hum i d l t  y 

50 

tt  specimens 

1 

Entry  dune  ns 

YES 

Thresh  delay 

.c?24tM 

Ex  tons  limit 

.90405 

1  9 . 4870 

D: 

.  0394 

r  :  MU 

Test  end 
Max  load 

Max  ox  t e 

status  : 
point  It 
ns  pnt  *1 

1" 

314 

6  4  '♦ 

1  Q .  /-noy. 

•*'.>  Ds 

.  o:tv 

37V 1  f  :  I 

Au::  t  1  1  np  ut  A r  r  a  y  11  1 

♦  *  *  *  +  *  *  «**#««»  »*•*««*»  ««**-«** 

Aux  i  1  i  ary  Input  Array  It  3 

**-»*«««  «-*«««  4  Y  «.»#««««  «*>«»«» 


»  *  «  *  4  #  *  *  *  ¥  *  *  4  «  »  *  «  «  4  *  «»««»»»  »»»«»«« 

»««*«+  «  M*  »««  '  «  *  4  •  m  *  * 
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CONFOSITES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
NONTEfiEY  CA  93990 


GRAPHITE  BUNDLE  TEST  (50  CM) 


Test  type:  Yern/F iber 

Operator  nane:  Nark  Jones 

Saaple  Identification:  100906 
Interface  Type:  9200  Series 
NacOin*  Pa'aeeters  of  test: 

Saeple  Rate  (pts/secl:  6.67 
Crosshead  Speed  (««/«in  ):  5.000 


Instron  Corporation 

Series  II  Automated  Naterials  Testi 

Test  Date:  Septeit.r  iv,i9B3 

Saaple  Type:  AS-L 


Huaidity  I  J  1:  50 

Teaperature  Ideg.  FI:  73 


'BATCH  Id'  100706 

’OIL  OR  NO  OIL1  OIL 


Cieer.sions: 


Spec.  1 


Lm.  Density  (te«>  £036 . 0 

Gauge  length  I  nr)  500.00 

Gut  t-f  1  spec  leers,  0  eicluded. 


i  Strain 

Load 

at 

at 

at 

Specimen 

Haxiiu« 

Kaurje 

K&xiflui 

Ntinbpr 

(n) 

txi 

UN) 

i 

3.640 

.7320 

.  J3c? 

Pean: 

3.660 

.7320 

.1329 

Standard 

Deviation: 


Keen  - 

2.00  •  Sd»:  --- . - . . 

Near  ♦ 

2. CO  *  Sdv:  . - . - . . 

Hiniaua:  3.660  .7320  .1327 


Nanrum:  3.660 


.7320 


.1327 


COMPOSITES  LAPQRATORV 
NAVAL  POSTGRADUATE  SCHOOL 
BCKTEPEY  CA  ?3?40 


6RAFHITE  BUNDLE  TEST  135.0  CHI 


Test  type:  Yarn/Fiber 

Operator  name:  Hark  Jones 


Instron  Corporation 

Series  II  Automated  Materials  Testing  5ystei  v4.01c 
Test  Cate:  September  10,1909 


Sa#ple  Identification:  10090? 

Interface  Type:  4300  Senes 
Machine  Parameters  of  test: 

Sample  Pate  (pts/sec):  5.00 
Crosshead  Speed  (RR/*in  1:  3.500 


Sasple  Type:  AS-4 

Humidity  I  X  ):  50 

Temperature  (deg.  f I :  73 


’PATCH  ID'  1009C7 

'OIL  OR  NO  OIL’  Oil 


Dieensions: 

Spec.  1 

Lin.  Density  <t?x)  3036.0 

Gauge  length  (sal  350.00 

Out  I  specimens,  0  excluded. 


Oisplceent  X  Strain  Load 
at  at  at 

Specimen  Naxiiut  Haxirua  Ha«»tu» 

Hvit5er  tt»l  \X)  t>N) 


I  3.370  .9480  .1710 


Mean  3.370  .94=0  .  1710 

Standard 

Dev  vat  ion:  . —  -  - 


Mean  - 
3.00  »  Sdv: 

Ht-'an  ♦ 

3.00  *  Sdv: 


Hinirje.: 

3.370 

.99  BO 

.1718 

HaxiiuR: 

3.370 

,9930 

.1718 
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VERSION  4.01c 


Sample  id 

1 00907 

Test  date 

September  10,1900 

Version 

A  .  0 1  c 

Version  date 

06/09/08 

Report  filett 

15 

Operator 

Mart  Jones 

X  conversion 

. 03937000 

X 

A/D  effse t 

.  00 '.“  i 

Y  conversion 

5.50465300 

Y 

A/D  offset 

.  0000 

Sample  rate 

5 . 00 

E:i 

tensometer 

NO 

A/D  range 

0 

Autostar  t 

OFF 

Cai lb  type 

AUTOMATIC 

Geomotr y 

CYLINDRICAL 

Ca  1  i  b  1 oad 

110.5311 

Ca  lib  e>:  t  ens 

.  0000 

Temperature 

73 

Hum i d i t  y 

50 

Test  type 

liTNICR 

tt 

spec i mens 

1 

Dnr  type 

AS  — 4 

Entry  d  linens 

YES 

Prfeil'  cherj 

1  0 .  00<’*  i*»0 

T hresh  delay 

.55401 

Load  limit 

1 15.40450 

E.“ 

tens  limit 

.  47/2  1  n 

Sample  cj  i  ii'rnn  i  cor.  i 

A  :  . 0394 

}< :  3f3324.oO0'> 

C: 

9. 

0453  III 

.0394  r :  NH 

SIpecimen 

1 

Test  end  status  :  10 

M<v:inium  load 

."IB .  573 

Ma 

load  p  (.'in 

i\  :  r-.r-j 

M,i::  r:::trnr. 

.  1  79 

Ma 

::  r ::  t  ci  ir.  pnt  4*  :  547 

Nijmbr*r  of  elr-innnbs  :  547 

Specimen  dimensions  : 

A:  .0093701  I« :  10324  .  OOOuoOO  C: 


Am ::  i  J  i  a r  y  1  npu  t  A r  r  a y  tt  1 

******  *  ♦  *  «  4  «-  #  #  *<*•**■**  *  «  **  V  *  *  ♦*■«***•»  * 

At  i  1  \  ary  Input  Array  t*  2 


9.0425200  D 


0393701  Ms  NO 


I 


■ 


* 

■ 


COMPOS: TES  LABORATORY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93940 


GRAPHITE  BUNDLE  UST  (25.0  CM) 


Test  type:  Yarn/Fiber 

Operator  name:  joe 


Instron  Corporation 

Series  fl  Automated  Materials  Testing  System  v4.0!t 
Test  Date:  Septeiber  1O,190B 


Sample  Type:  AS-4 


Humidity  (  X  ):  50 

Temperature  (deg.  f):  73 


'BATCH  ID’  100903 

'OIL  OR  NO  OIL’  NO  OIL 

Dimensions: 

Spec.  I 

Lin,  Density  (tex)  203b. 0 

Gauge  length  (ms»)  250.00 

Out  of  1  specirens,  0  excluded. 


Sample  Identification:  100908 
Interface  Type:  4200  Series 
Machine  Parameters  of  test: 

Sample  Rate  (pts/sec):  5.00 
Crosshead  Speed  (re/min  ):  2.500 


Specimen 

Nuaber 

Displcter.t 

at 

Maximum 

(mm) 

i  Strain 

at 

Maximum 

;xi 

Load 

at 

•Maximum 

(KN> 

1 

£.3’(i 

.9540 

.  1494 

fear. 

£.33'.' 

.9540 

.16-* 

Standard 

Mean  - 

2.00  *  Sdv: 

Mean  ♦ 

Minirjt: 

£.390 

.9540 

.1494 

Maxirjm: 

£.390 

.9540 

.1494 

r* 
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VERSION  4.01c 


Sample  id 

:  100900 

Test  date  : 

SeptcmbRi  1 0 . 190H 

Vf?r  si  on 

:  4.01c 

Version  date  : 

06/09/QO 

Report  filett 

:  15 

Operator  : 

joe 

X  conversion 

:  .03937000 

X  A/D  offset  : 

,00"0 

V  conversion 

:  c2  •  2 04  &2  3<  j  0 

Y  A/D  offset  : 

.  0*  M  Ml 

Sample  rate 

5 . 00 

E>:  tensometer  : 

NO 

A/D  range 

•  0 

Autostart  : 

OFF 

Ca 1  lb  type 

:  AUTOMATIC 

Geometry  : 

CYLINDRICAL 

Ca  1  ib  1 a Ad 

:  110.2311 

Ca lib  e:: tens  : 

.  0000 

T  empprature 

:  73 

Hum i d i t  y  : 

50 

Test  type 

:  DENIER 

#  specimens  : 

1 

Bar  type 

:  AS-4 

Entry  dimens  : 

YES 

Brea I.  chpcl 

•  10. 00000 

Thresh  delay  : 

. 22401 

Loaf]  limit 

:  112.40450 

Ex  tens  limit  : 

.49213 

Sariple  dimenr^ 

ions  : 

A :  .  0394 

f  i :  10324.0000 

C: 

9.0425  D: 

.  03V 4  F  :  Nil 

Specimen  \\ 

:  1 

Test  end  statu 

S  !  Ill 

Mm::  i  mum  1  (..-'ll 

:  :m."i  i 

1*1  a:;  load  point 

0  :  .  ■'  1 

l-la::  C'::  tens 

:  .174 

Ha::  ex  tens  pnt 

:  L.30 

N u  «nh  c?  r  r  •  f  e  1  i=-rr.  e  n  t  s  :  0  3 0 

'ippcinifi'ii  dimror-H-nr  : 

A:  .tX-iVSVvl  D :  1 8:).“*4  .  0«  M  #.;ii-iOO  C:  9.R4S^r«X'  D:  .0:v/37ol  L  :  MU 


Am::;  1  lfry  Input  Array  **  1 

*******  *«*-*-»-•»-  *******  *******  *******  *  *  *  «  *  *  ¥  *■**<•***  ******* 

Auxiliary  Input  Array  ft  c? 

«li  HMK  *  *  *  *  *  I  *  *******  *******  *******  #  *  *  *  *  4  *  *******  *******  .  }  M 
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